Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



METALLICITY CALIBRATIONS AND THE MASS-METALLICITY RELATION FOR STAR-FORMING 

GALAXIES 

Lisa J. Kewley^ 

Institute of Astronomy, University of Hawaii 



oo 

o 

O 



o 



> 

00 



Sara L. Ellison 

University of Victoria 
Draft version February 2, 2008 

ABSTRACT 

We investigate the effect of metalhcity cahbrations, AGN classification, and aperture covering frac- 
tion on the local mass-metallicity relation using 27,730 star-forming galaxies from the Sloan Digital 
Sky Survey (SDSS) Data Release 4. We analyse the SDSS mass-metallicity relation with 10 metalhc- 
ity calibrations, including theoretical and empirical methods. We show that the choice of metalhcity 
calibration has a significant effect on the shape and y-intercept(12-t-log(0/H)) of the mass-metallicity 
relation. The absolute metalhcity scale (y-intercept) varies up to A[log(0/II)] = 0.7 dex, depend- 
ing on the calibration used, and the change in shape is substantial. These results indicate that it 
is critical to use the same metalhcity calibration when comparing different luminosity-metallicity or 
mass-metallicity relations. We present new metalhcity conversions that allow metallicities that have 
been derived using different strong-line calibrations to be converted to the same base calibration. 
These conversions facilitate comparisons between different samples, particularly comparisons between 
galaxies at different redshifts for which different suites of emission-lines are available. Our new con- 
versions successfully remove the large 0.7 dex discrepancies between the metalhcity calibrations, and 
we reach agreement in the mass-metallicity relation to within 0.03 dex on average. We investigate 
the effect of AGN classification and aperture covering fraction on the mass-metallicity relation. We 
find that different AGN classification methods have negligible effect on the SDSS MZ-rclation. We 
compare the SDSS mass-metallicity relation with nuclear and global relations from the Nearby Field 
Galaxy Survey (NFGS). The turn over of the mass-metallicity relation at ~ 10^° M© depends on 
aperture covering fraction. We find that a lower redshift limit of z < 0.04 is insufficient for avoiding 
aperture effects in fiber spectra of the highest stellar mass (M* > 10^° M©) galaxies. 
Subject headings: galaxies: starburst — galaxies: abundances — galaxies: fundamental parameters — 
galaxies: spiral — techniques: spectroscopic 
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1. INTRODUCTION 

The relationship between metalhcity and stellar mass 
provides crucial insight into galaxy formation and evo- 
lution. Theory predicts that as time progresses, the 
mean stellar metalhcity of galaxies increases with age 
as galaxies undergo chemical enrichment, while the stel- 
lar mass of a galaxy will increase with tim e as galax- 
ies are built through merging processes f e.g.JPei fc Falll 
1995; Somerville & Primack 1999; Somer ville et al.ll2000l : 
Nagamine et al.ll2001t iCalura et al.ll20 04. and references 
therein). A correlation between mass and metalhc- 
ity arises if low mass galaxies have larger gas frac- 
tions than higher mas s galaxies, as is observed in lo- 
cal galaxies jMc Gaugh fc de Bloklll997l : iBell fc de Joni 
I2OO0I; iBoselh et al.ll200ll) . The detailed relationship be- 
tween metalhcity and mass may depend critically on 
galactic-scale o utflows driven by supernoy ae and stellar 
winds (see e.g.. lGarnettll2002t lPettinill2002l for a review). 
Thus, robust measurements of the mass-metallicity (MZ) 
relation may provide important clues into the impact of 
galactic-scale winds on the chemical history of galaxies. 

The MZ relation was first observed in i rregu- 
lar and blue compact galaxies (jLequeux et al.l 119791 : 
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iKinman &: DavidsonI |1981D . In subsequent work, lu- 
minosity was often used as a surrogate for mass be- 
cause obtain ing reliable mass es timates for galaxies was 
non-trivial. iRubin et al.l (|1984[ ) provided the first evi- 
dence that metalhcity is correlated with luminosity in 
disk galaxies. Further investigations solidified the cor- 
relation betwe en l uminosity and metalhcity in nearb\ 
disk galaxies dBothu n et al.l 11984 IWvse fc Silk 198E 



Skillman et al 



Zaritskv et al. 



19891 IVila-Costas fc Edmunds! 



199j;|Garnetdl200^. However, optical lu 



minosity may not be a reliable surrogate for the stel- 
lar mass of a galaxy because optical luminosities are 
sensitive to the level of current star formation and are 
extinguished by dust. Near infrared luminosities can 
be influenced by the age of the stellar population of 
a galaxy. Fortunately, reliable stellar mass estimates 
are now possible, thanks to new sta te-of-the-art stel 



lar evolutionary synthesis models (e.g..lSilva et al.l 



Leitherer et al 



•y syntnesis models le.g..lbiiva et al. 
1 Il999t [Fioc fc Rocca-Volmerangd 



1998 



1999; 



Bruzual fc Charlotll2003n" 

Key insight into the mass-metallicity relation has 
recently been obtained with large spectroscopic sur- 
veys such as the Sloan Digital Sky Survey (SDSS) 
and the 2 d egree Field Galaxy Redshift Survey (2dF- 
GRS) (e.g., iBaldrv et all 120021 : iSchulte-Ladbeck et all 



2003t iLamareille et all |2004 iTremonti efan 12004 
Gallazzi et al.ll2005D. Using the SDSS stellar masses, 



Tremonti et al.l (|2004l ; hereafter T04) characterized the 
local MZ relation for ~ 53,000 local galaxies. The MZ 
relation is steep for masses < 10^°'^ Mq and flattens 
at higher masses. T04 use chemical evolution models 
to interpret this flattening in terms of efficient galactic 
scale winds that remove metals from low mass galaxies 
(M < 10^°'^ M0). Hierarchical galaxy formation models 
that include chemical evolution and fccdbac k processes 
can reproduce the observed MZ relation (Dc Lucia et al.l 
l2004HDe Rossi et al.ll200llFinlator fc Davdl2007t ) . How- 
ever, these models rely on free parameters, such as feed- 
back efficiency, that are relatively unconstrained by ob- 
servations. Alternative scenarios proposed to explain 
the MZ relation include low star formation efficien- 



(jBrooks et al.ll2007 


) , and a variable intej 


tial mass function 


KoDDen et al.ll2007D 



The advent of large 8-10 m telescopes and effi- 
cient multi-object spectrographs enables the luminosity- 
metallicity (LZ) and, in some cases, the mass- 
metallicity relation to be char acterized to hig h red- 
shifts (|K obulnickv ~eraI1 ITOOl iCarollo fc Lillvl [200l^ 



Maier et al 



Pettini et al. 2001: Lillv et al.' 2003: Kobulnickv et all 
2003: Kobulnickv fc Kewley 2004; Shaolcv et al. 2004; 

Liang et al. 2004; Maier et al. 2005; 
Savaglio et al. 2005; Mouhcin e et al.. 
2006; "Maier et al. 2006; Lian g et all 
Evolution in the LZ and MZ relations 



2005 



Hovos et al 
20061 lErb et al 



2006fl . 



2004; 



are now predicted by semi-analytic models of galaxy 
formation within the A cold dark matter frame- 
work that include chemical hydrod ynamic simulations 
(iDe Lucia et al.ll2004l:lTisser a ct al. 2 005 iDe Rossi et al.l 
l2006HDave fc Oppenheimeil2007i ) . Therefore reliable ob- 
servational estimates of the LZ and MZ relations may 
provide important constraints on galaxy evolution the- 
ory. 

Reliable MZ relations require a robust metallicity cal- 
ibration. Common metallicity calibrations are based on 
metallicity-sensitive optical emission-line ratios. These 
calibrations include theoreti cal methods based on pho- 
toionization models (see e.g., iKewlev fc Dopital[2002i . for 
a review), empirical methods based on meas urements 
of th e electron-temperature of the gas, (e.g., iPilyuginI 
I2OOII: .Pcttini fc Pagcl 2004), or a combination of the 
two (e.g., IDenicolo et al.l 120021 ). Comparisons among 
the metallicities estimate d using these methods reveal 
large discrepancies ( e.g.. iPilvuginI 120011 : iBresolin et al.l 
[2OOI iGarnett et al.ll2004b[ ). These discrepancies usu- 
ally manifest as a systematic offset in metallicity esti- 
mates, with high values estimated by theoretical cali- 
brations and lower metallicities estimated by electron- 
temperature metallicities. Such offset s are found to be 
as large as 0.6 de x in log(0/H) units (|Liang et al.l 120061 : 
lYin et al.l[2007bh and may significantly affect the shape 
and zero-point of the mass-metallicity or luminosity- 
metallicity relations. 

Initial investigations into the extent o f these dis- 
crepan c ies have r ecently been r nade by [L iang et ajj 
(I2006D: lYi^ ejall (|2007bl lal). and iNagao et al. (2006^ 
iLiang et al.l ( 2006( 1 applied the metallicity calibrations 
from four authors to ~ 40, 000 galaxies from the SDSS. 
They showed that calibrations based on electron tem- 



perature metallicities produce discrepant metallicities 
when compar ed with ca l ibratio ns based on photoioniza- 
tion models. lYin et all (|2007b[ ) compared the theoreti- 
cal metal licities derived b y T 04 with Tg-based metallic i- 
ties from Pilyugin (2001) and 'Pilyugin & Thuanl 12003). 
They found a discrepancy of A[log(0/H)] = 0.2 dex be- 
tween the two Tg-based metallicities, and a larger dis- 
crepancy of A[log(0/H)] = 0.6 dex between the Tg-based 
methods and the theoretica l metho d. Similar results 
were obtain e d by lYin et al.l ()2007aD who extended the 
iLiang et al.l (|2006D work to low metallicities typical of 
low mass galaxies. 

The cause of the metallicity calibration discrepancies 
remains unclear. The discrepancy has been attributed 
to eithe r an unknown pr oblem with the photoionization 
models (jKennicutt et al. 2003) , or to temperature gradi- 
ents or fluctuations that may cause metallicities based on 
the electron tem perature ruethod to underest imate the 
true metalhcities ()Stasihskall2"00l l2005l : fBresol in 2006a). 
Until this discrepancy is resolved, the absolute metallic- 
ity scale is uncertain. 

The metallicity discrepancy issue highlights the need 
for an in-depth study into the effect of metallicity calibra- 
tion discrepancies and other effects on the MZ relation. 
In this paper, we investigate the robustness of the lo- 
cal MZ relation for star-forming galaxies. We focus on 
three important factors that may influence the shape, 
y-intercept, and scatter of the mass-metallicity relation: 
choice of metallicity calibration, aperture covering frac- 
tion, and AGN removal method. Our sample selection is 
described in Section [21 We describe the stellar mass esti- 
mates and metallicities in Sections [3l and [H respectively. 
We compare the mass-metallicity relation derived using 
10 different, popular metallicity calibrations in Section[5l 
We find larger discrepancies between the MZ relations 
derived with different metallicity calibrations than pre- 
viously found. We calibrate the discrepancies between 
the different calibrations using robust fits, and we pro- 
vide conversion relations for removing these discrepan- 
cies in Section [6l We show that our new conversions 
successfully remove the large metallicity discrepancies in 
the MZ relation. We investigate the effect of different 
schemes for AGN removal in Section [7l and we determine 
the effect of fiber covering fraction in Section [H We dis- 
cuss the impact of our results on the mass-metallicity and 
luminosity-metallicity relations in Section[9] Our conclu- 
sions are given in Section 1101 In the Appendix, we pro- 
vide detailed descriptions of the metallicity calibrations 
used in this study and worked examples of the application 
of our conversions. Throughout this paper, we adopt the 
flat A-dominated cosmology as me asured by the WM AP 
experiment (h = 0.72, = 0.29: ISpergel et al.|[20"03f) ). 

2. SAMPLE SELECTION 

We selected our sample from the SDSS Data Release 4 
(DR4) according to the following criteria: 

1. Signal-to-noise (S/N) ratio of at least 8 in 
the strong emission-lines [O II] AA3726, 9, 
H/3, [O III] A5007, Ha, [N II] A6584, and 
[S II] AA6717,31. A S/N> 8 is required for reli- 
able metallicity esti mates using estab l ishe d metal- 
licity cahbrations ([Kobulnickv et al.1 [l999i) . For 
each line, we define the S/N as the ratio of the 
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statistical error on the flux to the total flux, where 
the statistical errors are calculated by th e SDSS 
pipeline described in lTremonti et all (|2004D . 

2. Fiber covering fraction > 20% of the total photo- 
metric g'-band light. We use the raw DR4 fiber and 
Petrosian magnit udes to calcula te the fiber cover- 
ing fraction. Kewlev et al.l ()2005D found that a flux 
covering fraction > 20% is required for metallicities 
to begin to approximate global values. Lower cov- 
ering fractions can produce significant discrepan- 
cies between fixed-sized aperture and global metal- 
licity estimates. A covering fraction of > 20% 
corresponds to a lower redshift limit of z > 0.04 
for normal star-forming galaxies observed through 
the 3" SDSS fibers. Large, luminous star-forming 
galaxies larger redshifts to satisfy the covering frac- 
tion > 20% requiremend. We investigate residual 
aperture effects in Section [Sj 

3. Upper redshift limit z < 0.1. The SDSS star- 
forming sample becom es incomplete at r edshifts 
above z > 0.1 (see e.g.. lKewlev et al.ll2006D . With 
this upper redshift limit, the median redshift of our 
sample is z ~ 0.068. 

4. Stellar mass estimates m ust be available. Stellar 

Kauffmann et al.l (|2003[ ) 



mas ses were derived by 



and lTremonti et all (|2004l) 



We remove galaxies containing AGN from our sam- 
ple using th e opti cal classification criteria given in 
iKewlev et al.l ()2006( ). This classification scheme utilizes 
optical strong line ratios to segregate galaxies contain- 
ing AGN from galaxies dominated by star-formation. 
A total of 84% of the SDSS sample satisfying our se- 
lection criteria are star-forming according. This frac- 
tion of star-fo rming galaxies differs from the fraction in 
IKewlev et al.l pOOQ ) because we apply a more stringent 
S /N cut which removes many LINERs prior to classifica- 
tion. In Section [21 we investigate different AGN classifi- 
cation schemes and their effect on the shape of the MZ 
relation. 

The resulting sample contains 27,730 star-forming 
galaxies and does not include duplicates found in the 
original DR4 catal og. We note that our sample is smaller 
than the Trcmonti et al.l ([2004) sample because we ap- 
ply a more stringent S /N criterion and a stricter redshift 
range (T04 apply a redshift range of 0.005 < z < 0.25). 

We use the publically available emission-line fluxes 
that were calculated by the MPA/JHU group, (described 
in Tremonti et al. 2004). These emission- line fluxes were 
calculated using a sophistcated code that is optimized 
for use with the SDSS g alaxy spectra. This c ode a pplies 
a least-squares fit of the iBruzual fc CharlotI (|2003D stel- 
lar population synthesis models and dust attenuation to 
the stellar continuum. Once the continuum has been re- 
moved, the emission-line fluxes are fit with Gaussians, 
constraining the width and velocity separation of the 
Balmer lines together, and similarly for the forbidden 
lines. 

We correct the emission-line flu xes for extinct i on us - 
ing the Balmer decrement and the iCardelli et al.l (|1989f ) 
reddening curve. We assume an Ry = Av/E(B — V) = 
3.1 and an intrinsic Hq;/H/3 ratio of 2.85 (the Balmer 



decrement for case B recombination a t T= lO^K and 
He - 102 - lO^cm-^; IOsterbrockl[T98l . A total of 539 
(2%) of galaxies in our sample have Balmer decrements 
less than the theoretical value. A Balmer decrement less 
than the theoretical value can result from an intrinsically 
low reddening combined with errors in the stellar absorp- 
tion correction and/or errors in the line flux calibration 
and measurement. For the S/N of our data, the lowest 
E(B-V) measurable is 0.01. We therefore assign these 
539 galaxies an upper limit of E(B-V)< 0.01. 

To investigate aperture effects (Section [8]), we com- 
pare the SDSS mass-metallicity relation with the 
mass-metallicity relation deri ved from th e Ne arby 
Field Galaxy Survey (NFGS) (jJansen et al.l [2000 a b). 
IJansen et all |2000b) selected the NFGS o bjectively 
from the CfAl red shift survey ([Davis fc Peebl es 1983; 
iHuchra et al.lll983l ) to a pproximate th e local galaxy lu- 
minosity function (e.g.. (Marzke et al.|[l99l . The 198- 
galaxy NFGS sample contains the full range in Hubble 
type and absolute magnitude present in the CfAl galaxy 
s urvey. 

IJansen et al.l ([2000aD provide integrated and nuclear 
spectrophotometry for almost all galaxies in the NFGS 
sample. The covering fraction and metallicities of the 
nuclear and global spectra for the NFGS are described 
in Kewlev et al. (2005). The nuclear i?26 covering frac- 
tion ranges between 0.4 - 72%, with an average nuclear 
covering fraction of 10 ± 11%^. The covering fraction of 
the integrated (global) spectra is between 52-97% of the 
B-band light, with an average of is 82±7%. 

We apply the same extinction correction and AGN re- 
moval scheme to our NFGS supplementary sample as ap- 
phed to the SDSS sample. In the NFGS sample, 121/198 
galaxies can be classi fied using the i r nar row emission- 
lines according to the IKewlev et al.l ([2006D classification 
scheme. Of these, 106/121 (88%) are dominated by their 
star-formation. The NF GS integrated metallicities have 
been published bv , Kewlev et al.l ([2004f ) for several metal- 
licity calibrations. 



3. STELLAR MASS ESTIMATES 



The SDSS stellar masses were derived Oy 

[Tremonti et"an ([2004D and [Kauffmann et al.l ([20031 ) 
using a combination of z-band luminosities and Monte 
Carlo stellar population synthesis fits to the 4000A 
break and the stellar Balmer absorption line US a- The 
model fits to the 4000A break and IISa provide powerful 
constraints on the star formation history and metallicity 
of each galaxy, thus providing a more reliable indicator 
of mass than a ssuming a simple mass-to-ligh t ratio and a 
iKroupai ([2001[ ) Initial Mass Function (IMF). lDrorv et all 
(j2004l ) recently compared these spectroscopic masses for 
~ 17000 SDSS galaxies with (a) masses derived from 
population synthesis fits to the broadband SDSS and 
2MASS colors, and (b) masses calculated from SDSS 
velocity dispersions and effective radii. They concluded 
that the three methods for estimating mass agree to 
within ~ 0.2 dex over the 10® — 10^^ Mq range. 

An alte rnative method for es timating mass was pro- 
posed by [liir&deJoni ((IMil). Bell et al. used stel- 
lar population synthesis models to compute prescriptions 

^ The error quoted on the covering fraction is the standard error 
of the mean 



derived by 
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for converting optical cofors and pho t ometr y into stel- 
lar masses assuming a scaled iSalneteH ([Ml) IMF. This 
method is useful when near-IR colors are not available 
and spectral S/N is insufficient for reliable 4000A break 
and 116a measurements. We calculate masses for the 
NFGS galaxies by combining 2MASS J-band magnitudes 
with the B — R colors (R. Jansen, 2005, private com- 
munication). For all filters, we use 'total' magnitudes, 
i.e. the integrated light based on extrapolated radial 
surface brightness fits. We apply a search radius of 5 
arcsec in the 2MASS database, resulting in matches for 
85/106 star-forming galaxies. We calculate s tellar masses 
usmg the models of IBell fc de JontJ ('200i|), as param- 
eterised by [Rosenberg et al . (2005). For the compari- 
son between the SDSS and NFGS stellar masses (Sec- 
tion [5]), we assume a Salpeter IMF, and apply factors of 
1.82 and 1.43 to the SDSS (Kroupa) and NFGS (scaled 
Salpet e r) stellar m asses respectively iKaufFmann et al] 
((200l : lBeirfcdeJo ng(200lf). ^ 

Recently, iKannappan fc Gawiseii (|2007f ) compared 
stellar masses derived using stellar population synthesis 
fits to the NFGS spectra with masses derived using sev- 
eral methods, including the Bell et al. method. Kannap- 
pan & Gawiser find that the Bell et al. 2001 stellar mass 
prescription gives stellar masses that are 1.5x the pop- 
ulation synthesis approach (see Kannappan & Gawiser 
figure Ih). In the SDSS MZ relation, where stellar mass 
is in log space, a factor of ~ 1.5 would result in a shift of 
~ 0.17 dex. We consider this shift when comparing the 
NFGS and SDSS MZ relations (Section [3). 

4. METALLICITY ESTIMATES 

Metallicity calibrations have been developed over > 3 
decades from either theoretical models, empirical cali- 
brations, or a combination of the two. We apply 10 dif- 
ferent metallicity calibrations to the SDSS to investigate 
the impact of the metallicity calibration on the MZ re- 
lation. We divide the 10 calibrations into four classes; 
(1) direct, (2) empirical, (3) theoretical, and (4) calibra- 
tions that are a combination of empirical and theoretical 
methods. The empirical, theoretical and combined cal- 
ibrations all use ratios of strong emission- lines, and are 
often referred to collectively as " strong- line methods" to 
distinguish them from the "direct" method based on the 
weak 04363 auroral line. 

In this paper, we investigate the direct method, five 
theoretical calibrations, three empirical calibrations , and 
one " combined" calibration. We briefiy discuss each class 
of calibration below. The equations, assumptions, and 
detailed description of each method that we use are pro- 
vided in Appendix |^ and summarized in Table [1] 

4.1. Direct Metallicities 

The most direct method for determining metallicities 
is to measure the ratio of the 04363 auroral line to a 
lower excitation line such as [O III] A5007. This ratio 
provides an estimate of the electron temperature of the 
gas, assuming a classical H Il-region model. The electron 
temperature is then converted into a metallicity, after 
correcting for unseen stages of ionization. This method is 
sometimes referred to as the Ionization Correction Factor 
(IGF), or more commonly, the "direct" method, or the 
"Te" method. Determining metallicity from the auroral 
04363 line is subject to a number of caveats: 



1. The 04363 line is very weak, even in metal-poor 
environments, and cannot be observed in higher 
metallicity galaxie s without very sens itive, high 
S/N spectra (e..g.. iGarnett et al.|[200 4b). 

2. Temperature fluctuations or gradients within high 
metallicity H II regions may cause electron tem- 
perature metalliciti es to be underestimated by as 
much as ~ 0.4 dex (|Stasinskall200l l2005t iBresolinI 
[2006ah . In the presence of temperature fluctuations 
or gradients, [O III] is emitted predominantly in 
high temperature zones where O^"*" is present only 
in small amounts. In this scenario, the high elec- 
tron temperatures estimated from the 04363 line 
are not representative of the true electron temper- 
ature in the H II region, leading to syst ematically 
low metallicity estimates (see reviews by IStasihslal 
|2005; Bres ofin 2006al). 

3. The Te met hod may underestimate g lobal spectra 
of galaxies. iKobulnickv fc Zaritskvl (119991 ) found 
that for low metallicity galaxies, the Te method 
systematically underestimates the global oxygen 
abundance of ensembles of H II regions. 

High S/N ratio spectra can overcome the weakness 
of the 04363 line, and alternative auroral lines such 
as the [N II] A5755, [S III] A6312, and [O II] A7325 
lines are o bservable at higher metallicities than the 
4363 line (iKennicutt et al.ll2003t iBresolin et "all 12004 

funett et al.' ' 2004'bf r The theoretical investigation by 
tasiiiska 2005) predicts that these lines can provide ro- 
bust metallicities up to ~ solar (12 -I- log(0/II)= 8.7; 
AUende Prieto et al. 2001), but they may underestimate 
the abundance at metallicities above solar if temperature 
fluctuations or gradients exist in the nebula. 

4.2. Empirical Metallicity Calibrations 

Because 04363 is weak, empirical metallicity calibra- 
tions were developed by fitting the relationship between 
direct Te metallicities and strong-line ratios for H II 
regions. Typical calibrations a re based on the opt i- 
cal fine ratios [N II] A658 4/Ha (iPettini fc Pageil l2004f). 
([O III]/H/3)/([N II] /Ha), (|Pettini fc Pagell 12004 : here- 
after PP04), or the "Ra.y ra t fo (([O II] A3727 +\0 III] 
AA49 5 9,5007)/H/3 (iPilvugin '2001*; 'Pilvugin fc Thuad 
[200l iLiang et al.l 2007; Yin et al. 2007b)). PP04 
fit the observed relationships between [N II] /Ha, 
[O III]/H/3/[N II] /Ha and metalficity for a sample of 137 
H II regions. We refer to the Pettini fc Pagel methods as 
empirical because 97% of their sample has Te metallici- 
tieS; 

iPilyuginl (; hereafter POl 1200 If ) derived an empirical 
calibration for R23 based on Te-metallicities for a sample 
of H II regions. This ca libration has been updated by 
IPilvugin fc ThuanI (|2005L ; hereafter P05), using a larger 
sample of H II regions. 

We refer to strong-line metallicity calibrations that 
have been calibrated empirically from Te metallicities 
in H II regions as "empirical methods". In this paper, 
we apply the commonly-used empirical calibrations from 
POl (revised in P05), and PP04. These calibrations are 
described in detail in Appendix |Al and summarized in 
Table [TJ These empirical calibrations are subject to the 
same caveats as the Te-method described above. 
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4.3. Theoretical Metallicity calibrations 

The lack of electron temperature measurements at high 
metallicity led to the development of theoretical metallic- 
ity calibrations of strong-line ratios using photoionization 
models. These theoretical calibrations are commonly 
and confusingly referred to as " empirical methods" . The 
use of photoionization models to derive metallicity cal- 
ibrations is purely theoretical, and the use of the term 
"empirical" is a misnomer. We refer to photoionization 
model-based calibrations as "theoretical methods". We 
refer to all calibrations that are based on strong line ra- 
tios (i.e. including empirical and theoretical methods, 
but excluding the Tg method) as "strong- line" methods. 

Current state-of-the-ar t photoionization models 
such as MAPP I NGS (jSutherland fc Dopital fl99l 
Groves et al.l[200l [20061) and CLOUDY (iFerland eTaP 
1998f ) calculate the thermal balance at steps through 
a dusty spherical or plane parallel nebula. The ion- 
izing radiation field is usually derived from detailed 
stellar population syn thesis models such as Starburst99 
(jLeitherer et al.l[l999t ). The combination of population 
synthesis plus photoionization models allows one to 
predict the theoretical emission-line ratios produced at 
various input metallicities. 

Photoionization models overcome the temperature gra- 
dient problems that may affect Tg calibrations at high 
metallicities because photoionization models include de- 
tailed calculations of the temperature structure of the 
nebula. However, photoionization models have their own 
unique set of problems: 

1. Photoionization models are limited to spherical or 
plane parallel geometries. 

2. The depletion of metals out of the gas phase and 
onto dust grains is not well constrained observa- 
tionally 

3. The density distribution of dust and gas may be 
clumpy. This effect is not taken into account with 
current photoionization models. 

Because of these problems, discrepancies of up 
to 0.2 dex exist among the various strong-line 
calibrations based on photoionization models (e.g., 
[Kewlcy & Dopita 2002; Kobulnicky & Kcwlcy 2004, and 
references therein). Systematic errors introduced by 
modelling inac curacies are usually estimated to be ~ 
0.1-0.15 dex (|McGaughlll99ll:lKewlev fc Dopital[200l) . 
These error estimates are calculated by generating large 
grids of models that cover as many HII region scenar- 
ios as possible, including varying star formation histo- 
ries, stellar atmosphere models, electron densities, and 
geometries. Differences between the model assumptions 
and the true HII region ensemble that is observed in a 
galaxy spectrum are likely to be systematic, affecting all 
derived metallicities in a similar manner. 

Since systematic errors affect all of the direct, empiri- 
cal and theoretical methods for deriving metallicities in 
high metallicity (12 -I- log(0/H)> 8.6) environments, we 
do not know which method (if any) produces the true 
metallicity of an object. Fortunately, because the errors 
introduced are likely to be systematic, relative metallic- 
ities between galaxies are probably reliable, as long as 



the same metallicity calibration is used. We test this 
hypothesis in Section [9l 

Many theoretical calibrations have been developed 
to convert metallicity-sensitive emission-line ratios into 
metallicity estimates. Co mmonly used line ratios include 
[N II] A6584/[0 II] A3727 (jKewlev fc Dopita,.2002 . ; here- 
after KD02) and ( [ II] A3727 -H[0 111] AA4959,5007)/H/3 
(iPagel et al.l[l979l: lMcGaugh|[l99l IZaritskv et al.lfl99l 
Kobulnick y fc Kewlevll2004L : hereafter M91, Z94, and 
KK04 respectively) . In addition to the use of specific line 
ratios to derive metallicities, theoretical models can be 
used to simultaneously fit all observed optical emission- 
lin es to derive a metallic ity probability distribution, as 
in iTremonti et"an (|2004 ; hereafter T04). T04 esti- 
mated the metallicity for SDSS star-forming galaxies sta- 
tistically based on theoretical model fits to the strong 
emission-lines [O 11], H/3, [O III], Ha, [N II], [S 11]. 

In this paper, we apply the M91, Z94, KK04, KD02, 
and T04 theoretical calibrations, described in detail in 
the Appendix lAl Many empirical and theoretical metal- 
licity calibrations rely on the double- valued ([O II] A3727 
+ [0 III] AA4959,5007)/H^ line ratio, known as "R23". In 
Appendix IATT| we derive the [N II] /Ha and [N II] /[O 11] 
values that can be used to break the R23 degeneracy in 
a model-independent way. 

4.4. Combined Calibration 

Some metallicity calibrations are based on fits to the 
relationship between strong-line ratios and H II region 
metallicities, where the H II region metallicities are de- 
rived from a combination of theoretical, e mpirical and/or 
the di rect Te method. For example, the iDenicolo et al.l 
(|2002l ; hereafter D02) calibration is based on a fit to 
the relationship between the Te metallicities and the 
[N Il]/Ha line ratio for ^ 155 H II regions. Of these 
155 H 11 regions, ^ 100 have metallicities derived us- 
ing the Te method, and 55 H II have metallicities esti- 
mated using either the theoretic al M91 R2.3 method, or 
an em pirical method proposed bv lDiaz fc Perez- Monterg 
(|2000t ) method based on the sulfur lines. We refer to cal- 
ibrations that are based on a combination of methods 
as "combined" calibrations. In this paper, we apply the 
D02 combined calibration (described in Appendix IX)) . 

5. THE MZ RELATION: METALLICITY CALIBRATIONS 

In Figure [T] we show the mass-metallicity relation 
obtained using each of the 10 metallicity calibra- 
tions. There are insufficient galaxies in the SDSS with 
[O III] A4363 detections to determine an MZ relation 
using the Te metallicities. For the strong-line methods 
(i.e. all methods except the direct Te method), the red 
line shows the robust best-fitting 3rd-order polynomial 
to the data. The blue circles give the median metal- 
licity within masses of log(M) — 0.2 M0, centered at 
log(M) = 8. 6, 8. 8,..., 11 Mq. Both methods of charac- 
terizing the shape of the MZ relations produce similar 
results. The parameters of the best-fit polynomials and 
the rms residuals of the fit are given in Table [21 

The different strong-line calibrations produce MZ re- 
lations with different shapes, y-axis offsets, and scatter. 
T04 interpret the flattening in the MZ relation above stel- 
lar masses log(M) > 10.5 Mq in terms of efficient galac- 
tic scale winds that remove metals from the galaxies with 
masses below log(M) < 10.5 Mq. A similar flattening is 



6 



observed for the majority of the theoretical techniques. 
However, the MZ relations calculated using metallicity 
calibrations based on [N II] /Ha (D02 and PP04 N2) 
flatten at lower stellar masses log(M) ^ 10 because the 
[N II] /Ha line ratio becomes insensitive to metallicities 
for log([N II]/Ha)> -1 (or 12 + log(0/H)> 8.8 in the 
D02 or PP04 [N II]/Ha-based metallicity scale). The 
[N II] /Ha calibrations cannot give metallicity estimates 
above 12 + log(0/H)> 8.8, even if the true metallicity is 
higher than 12 + log(0/H)> 8.8^ 

The P05 empirical method (jPilvugin fc ThuanI [20051) 
is relatively flat for all stellar masses; between 8.5 < 
log(M/M©) < 11, the metallicity rises only ^ 0.2 dex on 
average. The majority of the H II regions used by P05 
have Te metallicities that are based on the [O III] A4363 
line. Because the [O III] A4363 line may be insensitive 
to (or saturate at) a metallicity 12 + log(0/H)~ 8.6, 
the P05 calibration may give a weak MZ relation for the 
SDSS. Interestingly, the original POl calibration (green 
line in panel (9) of Figure [T]) gives a steeper MZ rela- 
tion than the updated calibration (P05; red and blue 
lines). The updated P05 relation also produces lower ab- 
solute metallicities by ~ 0.2 dex co mpared with the or ig- 
inal POl method, as pointed out bv lYin et al.l (|2007bf ) in 
their comparison between POl, P05, and T04 metallici- 
ties. This change may be caused by the different H II- 
region abundance sets that were used to calibrate the 
original POl method and the updated version in P05. 

The direct Tg method is available for only 546/27,730 
(2%) of the galaxies in our SDSS sample. The 
[O III] A4363 line is weak and is usually only observed 
in metal-poor galaxies. The SDSS catalog contains 
very few metal-poor galaxies beca use they are intrinsi- 
cally rare, compact and faint (e.g..lTerlevich et al.l 119911 : 
iMasegosa e£al] 11994 Ivan Zed [20001). Panel 10 of Fig- 
ure in shows that a total of 477 Tg metallicities is in- 
sufficient to obtain a clear MZ relation. Because we are 
unable to fit an MZ relation using Tg metallicities, we 
do not consider the Tg method further in this work. 

The scatter in the MZ relation is large for all metallic- 
ity calibrations; the rms residual about the line of best- 
fit is 0.08 - 0.13. The cause of the scatter in the MZ 
relation is unknown. Our comparison between the dif- 
ferent metallicity calibrations shows that differing ion- 
ization parameter among galaxies does not cause or con- 
tribute to the scatter. The ionization parameter is explic- 
itly calculated and taken into account in some metallicity 
diagnostics (KD02, KK04, M91), but we do not see a re- 
duction in scatter for these methods. A full investigation 
into the scatter in the MZ relation will be presented in 
Ellison et al. (in prep). 

We directly compare the best-fit MZ curves for the 
9 strong-line calibrations in Figure [2 including both 
POl and P05. The top panel shows the rms scatter 
in metallicity about the mean in mass bins of width 
Alog(M/MQ) ~ 0.2. The major difference between 
the MZ curves is their position along the y-axis. The 
curves with the largest y-intercept are all photoioniza- 
tion model based (KK04, Z94, KD02, T04, M91). Among 
these photoionization model metallicities, the agreement 
is ~ 0.2 dex. This agreement is within the margin of error 
typically cited for these calibrations 0.1 — 0.15 dex for 
each calibration). Some calibrations consistently agree 
to within 0.1 dex (e.g., KK04 and Z94; KD02 and M91). 
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Fig. 1. — The mass-metallicity relation using the 10 differ- 
ent metaUicity calibrations listed in Table [l] The red line shows 
the robust best-fitting 3rd-order polynomial to the data. The 
blue circles give the median metallicity within stellar mass bins 
of Alog(M/M0) = 0.2, centered at log(M/M, 7)) = 8.6, 8.8, .... 11. 
We us e the updated calibration of P05 given bv lPilvuein fc ThuanI 
II2005I ) in panel 9. The original POl calibration is shown as a solid 
green line in panel 9. 

Comparisons between metallicities calculated using these 
consistent methods, such as KD02 and M91, are likely to 
be reliable to within 0.1 dex. However, comparisons be- 
tween methods that show large disagreement (such as 
KK04 and P05) will be contaminated by the large sys- 
tematic discrepancy between the calibrations. 

The lowest curves in Figure [2| are the MZ relations de- 
rived using the empirical methods (i.e. POl, P05, and 
the two PP04 methods). These empirical methods are 
calibrated predominantly via fits of the relationship be- 
tween strong-line ratios and H II region Te metallicities. 
There is considerable variation among the y-intercept 
of these Tg-based MZ relations; the P05 method gives 
metallicities that are ~ 0.4 dex below the PP04 meth- 
ods at the highest masses, despite the fact that both 
methods are predominantly based on H II regions with 
Te-metallicities. At the lowest stellar masses, this dif- 
ference disappears. The difference between the empirical 
methods may be attributed to the different H Il-rcgion 
samples used to derive the calibrations. At the highest 
metallicities, the PP04 methods utilize four H Il-regions 
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Fig. 2. — The robust best-fit mass-metallicity relations calculated 
using the different metallicity calibrations listed in Table [T] except 
the Te-method. The top panel shows the rms scatter in metallic- 
ity about the best-fit relation for each calibration in 0.1 dex bins 
of stellar mass. The y-axis offset, shape, and scatter of the MZ 
relation differs substantially, depending on which metallicity cali- 
bration is used. 

with detailed theoretical metallicities. These detailed 
theoretical metallicities may overcome the saturation at 
12 + log(0/H)~ 8.6 suffered by [O III] A4363 Te metal- 
licities. The P05 calibration includes some H II regions 
with metallicities estimat ed with the alt ernative auro- 
ral [N II] A5755 line from iKennicutt et a l. (2003). The 
inclusion of these [N II] A5755 metallicities may over- 
come the |0 II I] A4363 saturation problem. However, 
IStasinskal pOOSD suggest that the use of the [N II] A5755 
line in dusty nebulae will still cause Tg metallicities to be 
underestimated when the true metallicity is above solar. 
Our SDSS sample has a mean extinction of E(B-V)~ 0.3, 
or Ay ~ 1. The extinction is a strong function of stel- 
lar mass; for the largest stellar masses (M> IQ-'^^'^M©), 
the mean extinction is large E(B-V)~ 0.5, or Ay ^ 1.6. 
Clearly, dust is important in SDSS galaxies, particularly 
at the highest stellar masses where the largest discrepan- 
cies exist between the theoretical methods and the P05 
Tg-methods. 

In addition to the large difference in y-intercept be- 
tween the different metallicity calibrations, Figure [5] 
shows that the slope and turn-over of the MZ relation 
depend on which calibration is used. Therefore, it is 
essential to compare MZ relations that have been cal- 



culated using the same metallicity calibration. In the 
following section, we derive conversions that can be used 
to convert metallicities from one calibration into another. 

6. METALLICITY CALIBRATION CONVERSIONS 

Comparisons between MZ relations for galaxies in dif- 
ferent redshift ranges are non-trivial. Different suites of 
emission-lines are available at different redshifts, neces- 
sitating the use of different metallicity calibrations. Be- 
cause of the strong discrepancy in absolute metallicities 
between different calibrations, the application of differ- 
ent calibrations for galaxies at different redshifts may 
mimic or hide evolution in the MZ relation with redshift, 
depending on which calibrations are used. Because the 
metallicity discrepancies are systematic, we can fit the re- 
lationship between the different metallicity calibrations 
in order to remove the systematic discrepancies and ob- 
tain comparable metallicity measurements for different 
redshift intervals. 

We calculate conversion relations between the strong- 
line metallicity calibrations by plotting each calibra- 
tion against the remaining 8 calibrations and fitting 
the resulting metallicity-metallicity distribution with a 
robust polynomial fit. We refer to these metallicity- 
metallicity plots as Z-Z plots. Rows 1-3 of Figure [3] give 
six representative examples of SDSS Z-Z plots for the 
strong-line calibrations. The Z-Z plots between all nine 
strong- line calibrations for various S/N cuts are available 
at http://www.ifa.hawaii.edu/^kewley/Metallicity . The 
blue dashed 1:1 line shows where the Z-Z distribution 
would lie if the two calibrations agree. The robust best- 
fit polynomial is shown in red, and pr gives the robust 
equivalent to the standard deviation of the fit. Small 
values of pr indicate a reliable fit to the data. 

The majority of the Z-Z relations are close to linear and 
are easily fit by a 1st, 2nd or 3rd order robust polyno- 
mial. However, the P05 calibration produces a very non- 
linear relation with a large scatter when plotted against 
all other metallicity calibrations. These non-linear re- 
lations are not easily fit even with a robust 3rd order 
polynomial and we cannot provide conversions that will 
reliably convert to/from the P05 method. For compar- 
ison, the bottom row of Figure [3] shows the same plots 
calculated with the original POl calibration. Although 
the scatter is less severe in these plots, the relations be- 
tween POl and other diagnostics remain non- linear and 
are not easily fit with a robust 3rd order polynomial. 

For all other diagnostics, the metallicities or metallic- 
ity relations can be converted into any other calibration 
scheme, using 



y = a + bx 



(1) 



where y is the "base" or final metallicity in 12+log(0/H) 
units, a — d are the 3rd order robust fit coefficients given 
in Table [31 and x is the original metallicity to be con- 
verted (in 12 -I- log(0/H) units). For Z-Z relations where 
a 2nd order polynomial produces a lower pr than a 3rd 
order polynomial fit, d is zero. 

The conversion coefficients given in Table [3] are based 
on the fit order that produces the lowest pr value in our 
sample. Some R23 calibrations require two fits; one 2nd 
or 3rd order fit for the upper R23 branch and one linear 
fit for the lower branch. In these cases, the coefficents of 
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the upper and lower branch fits are listed in Table [3] as 
left and right columns, respectively. 

In Table [31 we give the range in x over which our cal- 
ibrations are valid. Our polynomial fits are only tested 
within these ranges and may not be suitable for con- 
verting lower or higher metallicities into another scheme 
outside these limits. We provide worked examples for the 
use of our conversions in Appendix [Bl 

Figure |3] shows the application of our strong-line con- 
versions to the best- fit MZ relations in Figure [U exclud- 
ing P05. The calibration shown for each panel represents 
the "base" (final) calibration into which all other MZ 
curves have been converted. The remaining discrepancy 
between the converted MZ relation and the base MZ re- 
lation is an indicator of both the scatter in our Z-Z plots 
and how well the Z-Z relations are fit by a robust polyno- 
mial. In Table 51 we give the mean residual discrepancy 
between the converted MZ relations and the base MZ 
relation. 

Our conversions reach agreement between the MZ re- 
lations to within 0.03 dex on average. The most reli- 
able base calibrations are those with the smallest residual 
discrepancies. The residual discrepancies differ because 
some Z-Z relations have less scatter and/or are more eas- 
ily fit by a simple polynomial. The KK04, M91, PP04 
03N2, and KD02 methods have the smallest residual dis- 
crepancies and are therefore the most reliable base cali- 
brations to convert other metallicities into. 

7. THE MZ RELATION: AGN REMOVAL METHODS 

The nebular emission line spectrum is sensitive to the 
hardness of the ionizing EUV radiation. Metallicities 
calculated from spectra that contain a significant con- 
tribution from an AGN may be spurious because the 
commonly-used metallicity calibrations are based on the 
assumption of a stellar ionizing radiation field. The 
standard optical diagn ostic diagrams for cl assification 
were first proposed bv I Baldwin et al.l (|1981h . based on 
the fine ratios [N II] /Ha vs [O III]/H/3, [S II] /Ha vs 
[O III]/H/3, and [O I] /Ha vs [O III]/H/3. This classifica- 
tion sclienie_wasj;evised_W (1985|) 
and lVeilleux fc Osterbro"cM (|1987l : hereafter VPS?) who 
used a combination of AGN and starburst samples with 
photoionization models to derive a classification line on 
the diagnostic diagrams to^separate A GN fro m starburst 
galaxies. Subsequently, Kewlcy ct al] (|2001l ; hereafter 
KeOl) developed a purely theoretical "maximum star- 
burst line" line for AGN classification using the standard 
diagrams. This theoretical scheme provides an improve- 
ment on the previous semi-empirical classification by pro- 
ducing a more consistent classification line from diagram 
to diagram that significantly reduces the number of am- 
biguously classified galaxies. The "maximum starburst 
line" defines the maximum theoretical position on the di- 
agnostic diagrams that can be attained by pure star for- 
mation models. According to the KeOl models, galaxies 
lying above the maximum starburst line are dominated 
by AGN activity and objects lying below the line are 
dominated by star formation. 

Although objects lying below the maximum starburst 
line are likely to be dominated by star formation, they 
may contain a small contribution from an AGN. We cal- 
culate the maximum AGN contribution that would allow 
a galaxy to be classified as star-forming with the KeOl 
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Fig. 3. — Examples of the relationship between different metal- 
licity calibrations. The robust 2nd order polynomial of best fit is 
shown as a red solid line. The 1:1 line (blue dashes) shows where 
the metallicities would lie if the calibrations agree. The robust 
equivalent to the standard deviation of the fit (pr) are shown for 
each plot. This figure illustrates the typical variation in scatter and 
shape between different metallicity calibrations. Figures showing 
the relation s between all 9 strong-line metallicity calibration s are 
available at |htt p : / / www . if a. hawaii . edu/ ~kewley /Metallicity | 

line on all three standard diagnostic diagrams using the- 
oretical galaxy spectra. Our AGN model is based on 
the 12 -I- log(0/H) = 8.9 dusty radiat ion-pressure dom- 
inated models by (jGroves et al.ll200^ . We use a typ- 
ical AGN ionization parameter of \og{U) = —2 and a 
power-law index of a = —1.4. We investiga t e the suite 
of starburst n iodels from iKewley fc Dopital (|2002f ) and 
iDopita et al.l (|2000r ). The starburst model that allows 
the maximum contribution from an AGN while remain- 
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Fig. 4. — The robust best-fit mass-metallicity relations calcu- 
lated using the 8 difli'erent metallicity calibrations listed in Table [T] 
converted into the base metallicities shown using our conversions 
from Table |3] 

ing classified as star-forming is zero-age instantaneous 
burst model with ionization paramete r qr = 1 x 10'' cm/s 
and metallicity 12+1ok(0 /H)= 8.9 by (jKewley fc Dopital 
[200llDopita et al.ll2000l ). The AGN contribution in this 
model is ^ 15%. 

We use this model to calculate the effect of a 15% AGN 
contribution to the metallicity-sensitive emission-line ra- 
tios. The AGN model contributes substantially to the 
[O III]/H/3 line ratio but has only a minor effect on the 
[N II] /[O II] ratio. Therefore, the effect of an AGN con- 
tribution of 15% is small (< 0.04 de x) on me tallicities 
calcu lated using the [N II] /[O II] ratio (jKewley fc Dopital 
l2002f ). but larger (0.1 — 0.2 dex) on metallici ties calcu- 
lated with calibrations containing [O III] fe.g..lMcGaughl 
[l99lHZaritsk^ et al."1994': 'Kobu lnickv fc Kewlevll2004r ). 

Recently. iKewley et al. (2006, ; hereafter Ke06) de- 
fined a new classification scheme based on all three diag- 
nostic diagrams that separates pure HII region-like galax- 
ies from HII- AGN composites, Seyferts, and galaxies 
dominated by low ionization emission line regions (LIN- 
ERs). This new classifi cation scheme include s an empir- 
ical shift applied by K auffmann et al.l (|2003l ; hereafter 
Ka03) to the KcOl line for the [N II]/Ha vs [O III]/H/3 
diagnostic. This shift provides a more stringent removal 



of objects containing AGN, and we recommend its use 
for metaUicities calculated using R23. 

We investigate whether the AGN classification scheme 
affects the shape of the MZ relation in Figure [5l For 
each metallicity calibration, we show the MZ relation for 
the three classification schemes KeOl (black dotted line), 
Ke06 (red solid fine) and V087 (blue dashed fine). These 
three classification schemes define 89%, 84%, and 76% of 
our SDSS sample as star-forming, respectively. There 
is negligible difference (< 0.05 dex) among the SDSS 
MZ relations for the three classification schemes. We 
note that the contribution from an AGN may be more 
important for samples that contain a larger fraction of 
HII- AGN composite galaxies, or galaxies at high redshift 
where limited sets of emission-lines limit the methods for 
AGN removal. For these cases, we recommend the use 
of either the KD02 [N II] /[O II] metalficity calibration 
(useful for log([N II]/[0 II])> -1.2), the PP04 [N II]/Ha 
calibration, or the D02 [N II] /Ha calibration. None of 
these three calibrations depend on the AGN-sensitive 
[O III]/H/3 line ratio. 
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Fig. 5. — Comparison between the mass -metallicity re l ations 
for three different methods of AGN removal; IKewley et al.l 1)20061 ) 
(red solid line), [K ewley ct al. (2001) (black dotted line), and 
IVeilleux fc Osterbroc k (,1987) (V087; blue dashed hne). The 
method of AGN removal has little effect on the MZ relation, except 
at low stellar masses here the V087 method gives a flatter slope. 
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8. THE MZ RELATION; APERTURE EFFECTS 

Our SDSS sample was selected with g'-band covering 
fractions > 20% because this value is the minimum cov- 
ering fraction req uired for metalliciti es to approximate 
the global values (jKewlev et al.l 120051 ). A covering frac- 
tion of 20% corresponds to a median redshift of z ~ 0.04 
which is the lower redshift limit used by T04 for their 
SDSS MZ relation work. The median g'-band aperture 
covering fraction of our sample is only ~34%, although 
the range of g'-band covering fractions is 20-80% (Fig- 
ure [6]). 



steUar masses (M > IO^^Mq). The NFGS global MZ re- 
lation flattens at a metallicity that is ~ 0.1 — 0.15 dex 
smaller than the metallicity at which the SDSS relation 
flattens. This difference is not caused by metallicity cal- 
ibration errors because the difference in upper turn-off 
is observed with all strong-line metallicity calibrations. 
Furtherm ore, the difference of log(M) = 0.17 dex be- 
tween thel Bell fc de Jongj (|2001h stellar mass relation and 
the SDSS Bruzual & Chariot model stellar masses cannot 
account for the difference between the SDSS and global 
NFGS MZ relations. 
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Fig. 6. — The distribution of g'-band fiber covering fractions 
in our SDSS sample. Our sample was chosen to have covering 
fractions > 20%. The dotted line at the top of the figure indicate 
the median (34.2) covering fraction of our sample. 



Strong metallicity gradients exist in most massive late- 
type spirals; H II region metallicities decrease by an order 
of magnitude from the inner to the outer disk (see e.g., 
[Shields! [l99 J. for a review). These gradients may cause 
substantial differ ences betw een the nuclear and global 
metallicities. iKewlev et al.l ('2005) investigate the effects 
of a fixed size aperture on spectral properties for a large 
range of galaxy types and luminosity. They conclude 
that minimum covering fractions larger than 20% may 
be needed at high luminosities to avoid aperture effects. 
Therefore the SDSS MZ relation may be affected by the 
fixed size aperture at the highest luminosities or stellar 
masses. 

To investigate the effect of the small median SDSS 
covering fraction on the MZ relation, in Figure [71 we 
compare the SDSS MZ relation (red solid line) with the 
nuclear (black dot-dashed line) and global (blue dashed 
line) MZ relations of the Nearby Field Galaxy Survey 
(NFGS). We show the KD02 metalhcity calibration (left 
panel) and the M91 calibration (right panel) for all 
datasets. Similar results are obtained with the other 
strong-line methods. The SDSS MZ relation lies in- 
between the NFGS nuclear and global relations at high 
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Fig. 7. — (Top Panel) Comparison between the robust best- 
fit SDSS Mass-Metallicity relation (red solid line) and the best- 
fit relations to the Nearby Field Galaxy Survey (NFGS) nuclear 
(blue dashed line) and global (black dot-dashed line) relations. 
Metallicities were calculated using both KD02 (left panel) and 
M91 (right panel) calibrations, and stellar masses are given as- 
suming a Salpeter IMF. (Bottom Panel) Comparison between the 
SDSS Mass-Metallicity relation (red solid line) and the NFGS 
global (filled circle) and NFGS nuclear (unfilled circle) data. At 
high stellar masses (M* > IO^^Mq), the SDSS metallicities are 
~ 0.1 — 0.15 dex larger than NFGS global metallicities at the same 
stellar mass. 



The difference between the SDSS and NFGS nuclear 
and MZ relations is probably driven by two factors: (1) 
fixed-size aperture differences and (2) different surface 
brightness profiles. The NFGS nuclear sample has a 
smaller mean covering fraction than the SDSS sample 
(~ 10% c.f. ~ 34%), giving higher nuclear metallici- 
ties in the NFGS than for SDSS galaxies with the same 
stellar mass. In addition, the NFGS and SDSS samples 
have different surface brightness profiles (traced by their 
half-light radii). The NFGS sample has a slightly smaller 
mean half light radius than our SDSS sample 3.0 kpc 
c.f. ~ 3.4 kpc respectively). Ellison et al. (in prep.) 
show that for the SDSS, galaxies with small g'-band half 
light radii (i.e more concentrated emission) have higher 
metallicities at a given mass than galaxies with large half 
light radii (more diffuse emission). 

The difference in half-light radii between the SDSS and 
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NFGS samples can not explain the difference between 
SDSS galaxies and global NFGS MZ relations at high 
stellar masses > lO'^'^Mobecause (a) the larger mean 
half-light radii of the SDSS sample would bias the SDSS 
towards lower metalhcities than the NFGS (Ellison et 
al. in prep), and (b) the NFGS global aperture covering 
fraction (~ 82%) captures most of the NFGS B-band 
emission. The half-light radius is irrelevant when the 
spectrum captures 100% of the B-band light. 

We conclude that a g'-band covering fraction of ^ 20% 
(or lower redshift limit of z = 0.04) is insufficient for 
avoiding aperture bias in SDSS galaxies with stellar 
masses > lO-'^^M©. The mean covering fraction for 
M^, > lO^^M© galaxies is 30.6±0.1%. A larger mean cov- 
ering fraction is required to obtain a reliable MZ relation 
at M, > 101°Mq. 

9. DISCUSSION 

We have investigated the effect of metallicity calibra- 
tions, AGN removal schemes, and a fixed-size aperture 
on the MZ relation. The choice of metallicity calibra- 
tion has the strongest effect on the MZ relation. The 
choice of metallicity calibration changes the y-intercept 
of the MZ relation significantly; the discrepancy between 
the metallicity calibrations is as large as 0.7 dex at the 
highest stellar masses. This discrepancy corresponds to 
a difference of 0.5 to 2.6 x solar at the peak metallicity 
of our SDSS sample. 

The existence of a ~ 0.4 — 0.5 dex discrepancy be- 
tween the Te an d theoretical metallicities is well known 
dStasinskai 120021 : [Kcn nicutt et al. 2003; iGarnett et"an 
l2004bD . Our results show that the discrepancy is larger 
than previously thought. This discrepancy is often in- 
terpreted as an unknown problem with the photoioniza- 
tion codes used to calibrate the str ong line r atios with 
metallicitv (iKennicutt et al.. 2003: G arnett et al. 20043 
iTremonti et al.l l2004l ) . However, recent investigations 
indicate that the Te methods may underestimate the 
metallicity when temperature fluctuatio ns or gradients 
exist within the emission-line nebulae ( Stasihskal l2005t 
iBrcsolin 2006a). These fluctuations, and hence the effect 
on Te are expected be the strongest at the highest metal- 
licities. We conclude that for the metallicities spanned 
by the SDSS sample, it is not possible to know which (if 
any) metallicity calibration is correct. Until the metal- 
licity discrepancies are resolved, only relative metallicity 
comparisons should be made. 

Relative metallicity comparisons rely on the ability 
of strong-line calibrations to consistently reproduce the 
metallicity difference between any two galaxies. For 
example, if two galaxies have metallicities of 12 -t- 
log(0/H)= 8.4 and 12-|-log(0/H)= 8.9 using one metal- 
licity calibration, the difference in relative metallicities 
(0.5 dex) should be the same using any other metallicity 
calibration, even if the absolute metallicities differ from 
one calibration to another. We test how well the strong- 
line metallicity calibrations maintain relative metallici- 
ties by selecting 30, 000 random sets of two galaxies from 
our SDSS sample. We measure the relative metallicity 
difference between the two galaxies from each set for 
each metallicity calibration. We give the mean differ- 
ence in relative metallicity and rms residuals in Table [5l 
The mean difference in relative metallicity is < 0.07 dex 
for all strong-line metallicity calibrations. The rms scat- 



ter is < 0.15 dex for all calibrations. The P05 method 
gives more discrepant relative metallicities to the other 
strong-line methods, with relative metallicity differences 
between 0.08 - 0.14 dex rms (c.f. 0.02 - 0.11 dex rms). 
The best agreement between relative metallicities oc- 
curs between the three theoretical R23 calibrations (M91, 
Z94, KK04), with relative metallicities agreeing to within 
0.02 — 0.05 dex rms. The small difference and rms resid- 
uals in relative metallicities for all 9 strong-line calibra- 
tions indicates that comparisons between galaxy or H II- 
region metallicities can be reliably made to within ~ 0.15 
dex, as long as the same base metallicity calibration is 
used for galaxies or H II regions. Our metallicity conver- 
sions aid relative metallicity comparisons between differ- 
ent samples of galaxies at different redshifts by empiri- 
cally removing the discrepancy between each metallicity 
calibration. In practice, if relative metallicity differences 
between galaxies or between samples is < 0.15 dex, we 
recommend the use of two or more metallicity calibra- 
tions to verify that any difference observed is real, and 
not introduced by the metallicity calibration applied. 

The SDSS sample is insufficient for determining the 
cause (s) of the metallicity discrepancy problem. Several 
ongoing investigations into the metallicity discrepancies 
may help solve this problem in the near future. These 
investigations include tailored photoionization models, 
high S/N spectroscopy of luminous stars in the Milky 
Way and nearby galaxies, metal recombination lines, IR 
fi ne structure lines, and temperature ffuctuation studies. 

IGarnett et al.l ()2004al ) applied tailored photoionization 
models to optical and infrared spectra of the H II re- 
gion COM 10 in M51. They found that the COM 10 
metallicity derived from the electron temperature us- 
ing the infrared [O III] 88^m line agrees with the the- 
oretical metallicity computed with the latest vers ion of 
the C LOUDY v90.4 photoionization code (Fcrland et al.l 
Il998f l. This theoretical metallicity is a factor of 2 smaller 
than the metallicity calculated with the previous version 
of CLOUDY (v. 74) that uses older atomic data. Clearly 
the optical emission-line strengths in the photoionization 
models are very sensitive to the atomic data used. How- 
ever, this sensitivity cannot explain the discrepancy ob- 
served in Figure [2] because T04 used the same version 
of Cloudy as Garnett et al. In spite of the use of mod- 
ern photoionization models with more accurate atomic 
data, the T04 MZ relation lies significantly higher than 
the methods uti li zing T e metalhcities (P05,D02,PP04). 
iMathis fc Wood! (|2005f l used Monte Carlo photoioniza- 
tion models to show that different density distributions 
are not a significant sourc e of error on the th eoreti- 
cal abundances. Recently, lErcolano et al.l (j2007l ) used 
new 3D photoionization codes to investigate the effect 
of multiply non-centrally located stars on the tempera- 
ture and ionization structure of H II regions. Ercolano et 
al. suggest that the geometrical distribution of ionizing 
sources may affect the metallicities derived using theoret- 
ical methods. Further theoretical investigations into the 
model assumptions, as well as tailored photoionization 
model fits to multi- wavelength data of spatially resolved 
star-forming regions may yield clues into whether the 
theoretical models contribute to the metallicity discrep- 
ancy. 

High S/N spectroscopy with 8- 10m telescopes can now 
provide metallicities for luminous stars and planetary 
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nebulae in nearby galaxies th at can be compared with 
H II region metallicities (s e e | Przvbilla et alJ I2007L for 
a review). lUrbaneia et al.l (|2005f ) analysed the chem- 
ical composition of B-type supergiant stars in M33. 
They find that the supergiant metallicities agree with 
H II region abundances derived using the Tg- method. 
Simi lar results were ob tained for local dwarf galax- 
ies (jBresolin et al.l [20061 ). however other investigations 
require a correction for electron ten iperature fluctua- 
tions before agreement can be reached (|Sim6n-Diaz et all 
[2OOI . 

Metal recombination lines provide a promising in- 
dependent baseline for metallicity measurements be- 
cause metal recombin ation lines depend only weakly 
on T. (jBresolinll2006al see e.g.,). Metal recombination 
lines are weak, but they h ave been observed in H II 
regions in the Milky Way (lEsteban et al.l |2004 120051: 
iGarcia-Roias et al] 120051120061). and recently in nearby 
galaxies (|Esteban et al.ll2002l : iPeimbert et aLll2005D . Re- 
combination methods typically agree with theoretical 
methods (e.g., Brcs olin., 2006a) . and predict larger metal- 
licities (by 0.2-0.3 dex) than the Te method. When 
the Te metallicities are corrected for electron temper- 
ature fluctuations, agreement is reached between re- 
combinati on and Te m et hods (iPeimbcrt ct al. 20051; 
Garcfa-Roias et al.l [200^, 120061 : [Peimbcrt ct al. 200i; 
Lopez-Sanchez et al.ll2007l ). 

Measurements of electron temperature fluctuations in 
nearby H II regions can resolve the disagreement between 
strong-line theoretical methods, and electron tempera- 
ture methods (Garcia-Roi as et al. 2006: Brcsolin 2003) 
in most cases (see however iHagele et al.l l2006l) . More 
electron temperature measurements are needed to ver- 
ify these results, particularly for high metallicity H II 
regions where electron temperature fluctuation measure- 
ments are lacking. 

Despite these promising investigations, the metallic- 
ity discrepancy problem remains unsolved at the present 
time. Until the metallicity discrepancy problem is re- 
solved, absolute metallicity values should be used with 
caution. In Table [5l we show that the metallicity cal- 
ibrations maintain relative metallicities better than 
0.15 dex rms, with the majority of theoretical methods 
maintaining relative metallicities within 0.04 — 0.1 dex 
rms. Therefore, studies of relative metallicity differences, 
such as comparisons between different galaxy samples, or 
between individual galaxies or H II regions, can be reli- 
ably made. If the size of the differences observed between 
different samples or different objects is 0.15 dex or 
less, we recommend the use of at least two independent 
calibrations to verify that the difference is calibration- 
independent. The KD02 and PP04 methods give (a) 
low residual discrepancies in relative metallicities, and 
(b) low residual discrepancy after other metallicities have 
been converted into these two methods. For the metal- 
Hcity range of the SDSS sample, the KD02 and PP04 N2 
calibrations are also independent of small contributions 
from an AGN. Because the KD02 and PP04 03N2 meth- 
ods maintain robust relative metallicities and are good 
base calibrations, we recommend the use of these two 
methods when deriving relative metallicities. 

10. CONCLUSIONS 



We present a detailed investigation into the mass- 
metallicity relation for 27,730 star-forming galaxies from 
the Sloan Digital Sky Survey. We apply 10 different 
metallicity calibrations to our SDSS sample, including 
theoretical photoionization calibrations and empirical Te 
method calibrations. We investigate the effect of these 
metallicity calibrations on the shape and y-intercept of 
the mass-metallicity relation. Using 30,000 galaxy sets 
sampled randomly from our SDSS sample, we investigate 
how well the 9 strong-line calibrations maintain relative 
metallicities. We find that: 

• The choice of metallicity calibrations has the 
strongest effect on the MZ relation. The choice of 
metallicity calibration can change the y-intercept 
of the MZ relation by up to 0.7 dex. Until this 
metallicity discrepancy problem is resolved, abso- 
lute metallicities should be used with extreme cau- 
tion. 

• There is considerable variation in shape and y- 
intercept of the MZ relations derived using the em- 
pirical methods. We attribute this variation to the 
different H II region samples used to derive the em- 
pirical calibrations. 

• The relative difference in metallicities is maintained 
to an accuracy of 12 + log(0/H)~ 0.02 - 0.1 dex 
for 9/10 cahbrations, and to within 12-flog(0/H)~ 
0.15 dex for all 9 strong-line calibrations. For rel- 
ative metallicity studies where the difference be- 
tween targets or between samples is < 0.15 dex, we 
recommend the use of at least two different calibra- 
tions to check that any result is not caused by the 
metallicity calibrations. 

We use robust fits to the observed relationship between 
each metallicity calibration to derive new conversion re- 
lations for converting metallicities calculated using one 
calibration into metallicities of another, "base" calibra- 
tion. Wc show that these conversion relations success- 
fully remove the strong discrepancies observed in the MZ 
relation between the different calibrations. Agreement is 
reached to within 0.03 dex on average. 

We investigate the effect of AGN classification scheme 
and fixed-size aperture on the MZ relation. 

• AGN classification methods have a negligible ef- 
fect on metallicities derived using [N II] /[O II] or 
[N II] /Ha. AGN classification can affect metallici- 
ties derived with R23 by < 0.15 dex. For the SDSS 
sample, AGN classification methods make negligi- 
ble difference in the shape or y-intercept of the MZ 
relation. For samples containing a larger fraction 
of starburst-AGN composite galaxies, or samples 
where AGN removal is not possible, we recommend 
the use of [N II] /[O II] or [N II] /Ha metalHcity di- 
agnostics. 

• The median g'-band aperture covering fraction of 
our SDSS sample is 34.2%. This covering fraction 
is insufficient for metallicities to represent global 
values at high masses (M> W^^Mq). The Nearby 
Field Galaxy global MZ is 0.1 — 0.15 dex lower than 
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the SDSS MZ relation at M> IO^^Mq. There- 
fore, the metaUicity at which the SDSS MZ rela- 
tion turns over is dependent on both the choice of 
metallicity calibration, and on the aperture size. 

We recom mend that meta ll icities be converted into 
either the iPettini fc Pagell (|2004[ ) method or the 
iKewlev fc Dopital (120020 method to minimize any resid- 
ual discrepancies, and to maintain accurate relative 
metallicities compared to other calibrations. These two 
diagnostics have the added benefit that at high metallic- 
ities, the Kewley & Dopita [N II] /[O II] and Pettini & 
Pagel [N II] /Ha calibrations are relatively independent 
of the method used to remove AGN. 

Future work into tailored photoionization models, high 
S /N spectroscopy of luminous stars in the Milky Way and 



nearby galaxies, metal recombination lines, IR fine struc- 
ture lines, and temperature fluctuation studies may help 
resolve the metallicity discrepancy problem in the near 
future. Until then, only relative metallicity comparisons 
are reliable. 
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APPENDIX 



METALLICITY CALIBRATIONS: EQUATIONS AND METHOD 
Breaking the R23 Degeneracy 

Many empirical and theoretical metaUicity calibrations rely on the ([O II] A3727 -|-[0 III] AA4959,5007)/H/3 line 
ratio, known as "R23" • The major drawback to using R23 is that it is double- valued with metallicity; R23 gives both 
a low metallici ty estimate ("lower branch") and a high estimate ("upper branch") for most values of R23 (see e.g., 
iKobulnickv fc Kewley 2004 for a discussion). Additional line ratios, such as [N II]/Ha, or [N II]/[0 II], are required 
to break this degeneracy. 

The SDSS catalog contains very few metal-poor galaxies (jlzotov et al.l 120041 : iKniazev et al.l l2003l 120041 : 
iPapaderos et al.l [2006t llzotov et al.l l2006a[ ). Metal poor galaxies are of ten lacking in magnitude-limited emission- 
line surveys b ecause they are intrinsically rare, compact and faint (e.g.. 'Terle vich et al.lll99ll iMa scgosa et al."1994l: 
Ivan "Ze3 120001 ). For the purpose of investigating the uppe r and lower R 2 3 bra nches, we supplement the SDSS sample 
with (a) the low met allicity galaxy sample described in iKewlev et al.l (|2007[ ) and iBrown et al.l ()2006[ ). and (b) the 
iKong fc Chend (|2002| ) blue compact galaxy sample. 

Note that we do not calculate an initial metallicity from an [N II] /Ha or [N II] /[O II] metallicity calibration because 
in some cases, a systematic discrepancy between a metallicity calibration based on [N II] /Ha or [N II] /[O II] and the 
calibration based on R23 will cause galaxies to be improperly placed on the upper or lower branch of R23. For example, 
an [N II] /Ha metallicity calibration that systematically produces higher estimates than the subsequent R23 calibration 
may cause metallicities to be erroneously estimated from the upper R23 branch. 

We use the [N II] /[O II] ratio to break the R23 degeneracy for our SDSS sample. The [N II] /[O II] ratio is not sensitive 
to ionization parameter to within (±0.05 dex), and it is a strong function of metallicity above log([N II]/[0 II])> —1.2 
(jKewley & Dopita 2002). Figure [Hi shows that the division between the R23 upper and lower branches occurs at 
log([N II]/[0 II])~ —1.2 for the SDSS and supplementary samples. For comparison, Figure[5jD shows the theoretical re- 
lation ship between [N II] /[O II] and R23 using the population synthesis and photoionization models of lKewlev fc Dopital 

(|2002| ). The observed R23 peak at log([N II]/[0 II]) 1.2 corresponds to a metaUicity of 12-hlog(0/H)~ 8.4 according 

to the theoretical models. 

For galaxies at high redshift, the [N II] /[O II] ratio cannot be used to break the R23 degeneracy because either (a) 
[N II]/[0 II] cannot be corrected for extinction due to a lack of reliable Balmer line ratios, and/or (b) [N II] and [O II] 
are not observed simultaneously in a given spectrum. In this case the [N II] /Ha ratio is used (Figure [9]). Figure 
shows that the division between the R23 upper and lower branches occurs between —1.3 <log([N II]/Ha)< —1.1 for the 
SDSS and supplementary samples. The division between the upper and lower R23 branches using [N II] /Ha (Figure|9^) 
is less clear than for [N II] /[O II] (Figure EK) because the [N II] /Ha ratio is less sensitive to metallicity, and more 
sensitive to ionization parameter, than [N II] /[O II]. 

We check whether our empirical [N II]/Ha division between the upper and lower R23 branches —1.3 <log([N II]/Ha)< 
— 1.1 is compatable with our [N II] /[O II] division (log([N II] /[O II] )< —1-2) by comparing the number of galaxies placed 
on the upper and lower branches using each ratio. For log([N II]/Ha)< —1.3, the majority of galaxies (150/175; 86%) 
lie on the lower R23 branch according to their [N II]/[0 II] line ratios. For log([N II]/Ha)> —1.1, the upper R23 branch 
can be clearly seen. For —1.3 <log([N II]/Ha)< —1.1, [N II]/Ha cannot discriminate between the upper and lower 
R23 branches. Figure|9h) shows that galaxies with —1.3 <log([N II]/Ha)< —1.1 are likely to have metallicities that are 
close to the R23 maximum. For the SDSS galaxies with -1.3 <log([N II]/Ha)< -1.1 in Figure |9ji, the [N II]/[0 II] 
ratios indicate that 634/1127 (56%) of SDSS galaxies lie on the upper branch and 493/1127 (44%) SDSS galaxies lie 
on the lower branch. 

For the R23 methods in this paper, we use the [N II] /[O II] line ratio to break the R23 degeneracy. 
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Fig. 8. — (a) The observed relationship between the metaUicity-sensitive [N II] A6 584/[0 II] A3727 Hne ratio and the commonly-used 
([O II] A 3727 +[0 III ] AA49 59,5007)/H/3 ratio. The SDSS galaxies (black), the lKong & Cheng (2002) blue compact galaxy sample (blue) 
and the IBrown et al.l II2006I ) low metallicity galaxy samples (red) are shown. The [Nil] /[Oil] ratio is a strong monotonic function of 
metallicity to log([NII]/[OII])> —1.2, while R23 has a maximum at log(R23 )~ 0.9. For our samples, the R23 maximum is likely to occur at 
log([NII]/[OII])~ —1.2. This value can be used to break the R23 degeneracy for galaxies where [NII]/[OII] can be corrected for extinction 
using the Balmer Decrement, (b) The theoretical relationship between the [Nil] /[Oil] and R23 line ratios using the stellar population 
synthesis and photoionization model grids of Kewley & Dopita (2002). Models are shown for constant metallicities of 12 + log(0/H)= 
7.9,8.2,8.6,8.9, 9.1, 9.2 and ionization parameters of g = 1 X 10^, 2 X lO'', 4 X 10^, 8 X 10^, 1.5 X 10® cms'^. We choose a break between 
the R23 upper and lower branches at log([NII]/[OII])~ —1.2, which corresponds to a metallicity of 12 + log(0/H)~ 8.4 according to the 
theoretical models. 



Theoretical Photoionization Methods 
McGaugh (1991) - M91 

The iMcGaughl (Il991l) calibr ation of R23 is based on detailed H II region models using the photoionization code 
CLOUDY (|Ferland et al.l 119981 ). The M91 calibration includes a correction for ionization parameter variations. We 
use the [N II] /[O II] line ratio to break the R23 degenerac y, as described in Sectio n [A. 11 and we apply the analytic 
expressions for the M91 lower and upper branches given in iKobulnickv et al.l (|1999D : 

12 + \og{0/Yi)io^er = 12 - 4.944 + 0.767a; + Qm2x^ 

-j/(0.29 + 0.332a; - 0.331a;2) (Al) 



12 + \og{0/B)upper = 12 - 2.939 - 0.2a; - 0.237a;2 

-0.305a;3 - 0.0283a;'* - y(0.0047 
-0.0221a; - Q.lQ2x^ - 0.0817x^ 

-0.00717x'^) (A2) 

where X = log{R,,) = log ( fO"1^3727+[OIII]A4959+[OIII]A5007 ^) ^ ^ ^ ^^^(0^^) ^ |- [OIII]A4959+[OIII] A5007 ^ 

estimated accuracy of the M91 calibration is ^ 0.15 dex. 
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Fig. 9. — (a) The observed relationship between the metallicity-sensitive [N II] A6584/Ho hne ratio and the commonly-used ([O II] 
A37 27 +[0 III] AA4959 ,5007)/H/3 ratio. The SDSS galaxies (black), the JionglErChcnK (2002) blue compact galaxy sample (blue) and 
the [B rown et al.l I I2006I) low metallicity galaxy samples (red) are shown. For log([NII]/Ha) > —1.1, galaxies are likely to lie on the 
upper R23 branch. For —1.3 <log([NII]/Ha )< —1.1, [NII]/H« cannot discriminate between the upper and lower R23 branches. For 
log([N II]/Ha)< —1.3, galaxies lie on the lower R,23 branch, (b) The theoretical relationship between the [N II] /[O II] and R23 line 
ratios using the stellar population synthesis and photoionization model grids of Kewley & Dopita (2002). Models are shown for constant 
metallicities of 12 + log(0/H)= 7.9,8.2,8.6,8.9, 9.1,9.2 and ionization parameters of g = 1 X lO'', 2 X lO'^, 4 X lO'', 8 X lO'^, 1.5 X 10** cms~^. 



Kewley & Dopita (2002) - KD02 

The iKewlev fc Dopital (|2002l ) calibrations are based on a self-consistent combination of detailed stellar population 
synthesis and photoionization models. The estimated accuracy of the KD02 calibrations is ^ 0.1 dex. This estimate 
is derived by varying the major assumptions in the stellar evolution and photoionization models (including the star 
formation prescription, electron density, and the initial mass function) . KD02 outlined a "recommended" approach to 
deriving metallicities that uses the [N II] /[O II] line ratio for high metallicities and a combination of R23 calibrations 
for lower metallicities. We use a revised version of the KD02 calibration. For log([N II]/[0 II])> —1.2 , we use the 
original KD02 [N II] /[O II] metallicity calibration given by 

\og{[NII]/[OII]) = 1106.8660 - 532.15451Z + 96.373260^2 - 7.8106123^^ + 0.23928247^" (A3) 

where Z = \og{0 / H) + 12. We use the IDL task fz-roots to solve the 4th order polynomial for Z. The coefficients 
in Equation IA3I are based on the theoretical <? = 2 x 10*" cm/s relationship between [N II]/[0 II]and Z. However, 
the detailed relationship between [N II] /[O II] and Z is independent of ionization parameter to within ~ 0.1 dcx for 
log([N II]/[0 II])> —1.2 and the ionization parameters covered by the SDSS (g = 1 x lO'' — 8 x lO*" cm/s). 

For log([N II]/[0 II])< -1.2 (or KD02 12 -1- log(0/H)< 8.4), KD02 recommend using an average of R23 methods. In 
this regime, we use t he av erage of the KK04 lower banch R23 calibration (equation I A6p and the lower branch 1VI91 R23 
calibration (equation [Al]). Both of these calibrations correct for ionization parameter variations. 

Kobulmcky & Kewley (2004) - KKO4 

iKobiilnickv fc KewlevI (I2OOI 1 use the stellar evolution and photoionization grids from iKewlev fc Dopital (j2002f ) to 
produce an improved fit to the R23 calibration. The estimated accuracy of the KK04 method is ~ 0.15 dex. 

The R23 calibration is sensitive to the ionization state of the gas, particularly for low metallicities where the R23 
line ratio is not a strong function of metallicity. The ionization state of the gas is characterized by the ionization 
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parameter, defined as the number of hydrogen ionizing photons passing through a unit area per second, divided by 
the hydrogen density of the gas. The ionization parameter q has units of cm/s and can be thought of as the maximum 
velocity ionization front that a radiation field is able to drive through the nebula. The ionization parameter is typically 
derived using the [O III]/[0 II] line ratio. This ratio is sensitive to metallicity and therefore KK04 recommend an 
iterative approach to derive a consistent ionization parameter and metallicity solution. We first use the [N II] /[O II] 
ratio to determine whether each SDSS galaxy lies on the upper or lower R23 branch. We then calculate an initial 
ionization parameter by assuming a nominal lower branch (12 + log(0/H)= 8.2) or upper branch (12 + log(0/H)= 8.7) 
metallicity using equation (13) from KK04, i.e. 

log(g) = {32.81 - 1.153y2 

+ [12 + log(0/H)](-3.396 - 0.025y + 0.144^)} 
X {4.603- 0.3119?;- 0.163J/2 

+ [12 + log(O/H)](-0.48 + 0.0271?; + 0.02037?/2)}~i (A4) 

where y = log032 = log([OIII] A5007/[OII] A3727). The initial resulting ionization parameter is used to derive an 
initial metallicity estimate from KK04 equation (16) for log([N II]/[0 II])< -1.2 (12 + log(0/H)< 8.4), or KK04 
equation (17) for log([N II]/[0 II])> -1.2 (12 + log(0/H)> 8.4): 

12 + log(0/H);o^er = 9.40 + 4.65a; - 3.172;^ 

- log(g) (0.272 + 0.547a; - O.SlSa;^) (A5) 



12 + log(0/H)„pper = 9.72 - 0.777a; - QMlx^ - 0.072a;^ - O.Sllx^'a; 

- log(g) (0.0737 - 0.0713a; - O.Wla;^ + 0.0373a;3 - 0.058a;'*) (A6) 

where x = logRas = log >3727+[oin] aa4959,5007 ^ ^ Equations El and Ell (or El]) are iterated until 12 + log(0/H) 
converges. Three iterations are typically required to reach convergence. 

Zantsky et al. (1994) - Z94 

The IZaritskv et al.l ([1994') cal ibration is based on th e R23 line ratio. This ca li bration is der i ved fr om the average of 
three previous calibrations bv .Edmunds fc Pagell (|1984D ; iDopita fc Evans! (|1986f) ; iMcCall et "all (|1985f) . The uncertainty 
in the Z94 calibration is estimated by the difference in metallicity estimates between the three calibrations. Z94 
provide a polynomial fit to their calibration that is only valid for the upper R23 branch (i.e. 12 + log(0/H)> 8.4, or 
log([N II]/[0 II])> -1.2). 

12 + log(0/H) = 9.265 - 0.33x - 0.202a;2 - 0.207^^ - 0.333a;'* (A7) 

where x ~ logR23 — log ^^**^ A3727+[oin] aa4959,5007 ^ ^ solution for the ionization parameter is not explicitly 
included in the Z94 calibration. 

Tremonti et al. (2004) - TO4 

T04 estimated the metallicity for each galaxy statistically based on theoretical model fits to the strong emission- lines 
[O II] , H/3, [ O III], Ha, [N II], [S I I]. Th e model fits were calculated using a co mbination of stellar po pulation synthesis 
models from lBruzual fc CharlotI (|2003l ) and CLOUDY photoionization models H^and et all ()1998D . The T04 scheme 
is more sophisticated than the other theoretical methods because it takes advantage of all of the strong emission lines 
in the optical spectrum, allowing more constraints to be made on the model parameters. Calibrations of various line 
ratios to the theoretical T04 method are given by iNagao et al.l (120061) and Liang et al. (2006). We use the original 
T04 metallicities, available from http:/ / www. mpa-garching. mpg. de / SDSS/\ for this study. 

Te method 

We derive the gas-phase oxygen abundance followin g the proced ure outlined in llzotov et al.l (l2006bl). This pr o cedur e 
utilizes the electron-temperature (Tg) calibrations of lAUed (|1984f ) and the atomic data compile d bv IStasinskal j2005| ). 
Abundances are determined within the framework of the classical two-zone Hll-region model (jStasinskal 11980 ). The 
ratio of the auroral [O III] A4363 and [O III] AA4959, 5007 emission-lines gives an electron temperature in the 0++ 
zone. We derive electron densities measured using the [S II] A6717/[S II] A6731 line ratio. These electron temperatures 
are insensitive to small variations in electron density; we obtain the same Te with an electron density of 367 c m~^. 
The electron temperature of the O"*" zone is calculated assuming Tc(0+) = 0.7Te(O++) +0.3 (|Stasihskalll98l . We 
calculate the metallicity in the 0+ and 0++ zones assuming 



0/H = 0+/H+ + 0++/H+ 



(A8) 
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The uncertainty in the absolute 0/H metallicity determination by this Tg method is ~ 0.1 dex. This intrinsic 
uncertainty is the dominant error in our Tg metallicity determination, and i ncludes errors in the use of simplified H II 
region models and possible problems with electron temperature fiuctuations (jPagel fc Tautvaisienel[T997l ) . Fortunately, 
these errors affect all Tg-based methods in a similar way and the error in relative metallicities derived using the same 
method is likely to be << 0.1 dex. 

This "c lassical" T^ appr oach does not take unseen stages of ionization or electron temperature fluctuations into 
account. iBresolinI (|2006a[ ) notes that if electron temperature fluctuations are substantial and are not taken into 
account, Tg-based calibrations can only provide a lower limit to the true metallicity, particularly in the high metallicity 
regime where Tg fluctuations ap pear stronger. We find that the Tp method does not produce any SDSS metallicities 
of solar (12 + log(0/H)- 8.69; lAllende Prieto et all I2OOII : lAsplund et al.l 12004 ) or above, even for galaxies where 
the fiber only captures 20 — 30% of the central g'-band galaxy fight. Covering fractions o f 20 — 30% correspond to 
diameters of ~ 5 kpc for the mean size of nearby star- forming galaxies (jKewlev et al.l 120051 ). Spiral galaxies typically 
have metallicities within similar apertures that are ~ 1 — 2 x solar, measured using various independent methods (see 
iHenrv fc Worthe"vlll999l for a review). For example, our Gal axy has consistent centra l metallicities within the central 
5 kpc of 1 — 2 X solar from studies of planetary nebulae (jMartins fc Viegail200C ), IR fine stru cture lines in H II 



regions (jSimpson et al. 1995; Afflcrbach et al- 1997), and radio recombination lines ( Quireza et al.i i2006). 

We conclude that Tg metallicities should be used with caution when other line ratios (such as [N II] /Ha and 
[N II] /[O II]) indicate upper branch (Figures [8] and [9]) or supersolar metallicities. 

Empirical T,, fit methods 

H II regions with electron temperature-based metallicities have been used in many studies to derive empirical fits to 
strong-line ratios that can be applied to H II regions and galaxies where the 04363 line is not observed. 

Petttm & Pagel (2004 ) - PP04 03N2 & N2 

iPettini fc Pagell (j2004l ) derived two new methods for measuring metallicities in galaxies at high redshift. At high 
redshift, obtaining a reddening estimate is difficult and in some cases, impossible, and flux calibration in the infrared 
is non-trivial. Ratios of lines that are very close in wavelength do not require reddening correction or flux calibration. 
PP04 fit the observed relationships between [N II] /Ha, ([O III]/H/3)/([N II] /Ha) and the Tg-based metafiicity for a 
sample of 137 H II regions. Of these 137 H II regions, 131 have Tg-based metallicities and 6 high metallicity H II 
regions have metallicities derived using detailed photoionization models. Because the vast majority of H II regions in 
the PP04 sample have Tg-based metallicities, we refer to the PP04 method as an empirical Tg fit method. The fit to 
the relationship between Tg metallicities and the ([O III]/H/3)/([N II]/Ha) ratio is: 



12 + log(0/H) 8.73 - 0.32 x 03N2, (A9) 

where 03N2 is defined as 03N2 = log ( fJ!!]x658l/Ha ) ■ Equation El] is only valid for 03N2> 2. We refer to this 
calibration as "PP04 03N2". 

PP04 fit the relationship between Tg metallicities and the [N II] /Ha ratio by a line and a third-order polynomial. 
We use the polynomial fit given by 

12 -I- log(0/H) 9.37 + 2.03 x N2 + 1.26 x + 0.32 x N2^ (AlO) 

where N2 is defined as N2 = \og{[NII] X6584/Ha). Equation ElO] is valid for -2.5 < N2 < -0.3. 

Because the PP04 method was derived using a fit to H II regions rather than galaxies, we check whether the PP04 
relations are suitable for metallicity estimates of the SDSS sample. In Figure [10] we show the relationship between 
(a) N2 and Tg metallicities, and (b) 03N2 and Tg metallicities for the SDSS galaxies in our sample with measurable 
(S/N> 3) [O III] A4363 lines. The dashed line indicates the PP04 calibrations based on H II regions, while the dotted 
lines encompass 95% of the H II regions in the PP04 sample. The majority of the SDSS galaxies lie within the 
PP04 95 percentile fines. However, 47/546 (9%) and 69/546 (13%) of SDSS-Tg galaxies have Tg metalficities that 
lie below the 95 percentile line in the N2 and N203 diagrams, respectively. These galaxies have high [N II] /Ha and 
[N II]/[0 II] ratios (log([N II]/Ha)> -1; log([N II]/[0 II])> -0.8), indicating supersolar metalficities, according to aU 
of the [N II]/Ha and [N II]/[0 II]-based metalficity diagnostics. Both Figure[TO]and the high [N II]/Ha and [N II]/[0 II] 
ratios suggest that the Tg-method underestimates metallicities for galaxies that lie below the PP04 95 percentile line. 



Pilyugin (2005) - P0.5 

iPilvuginl ()2001l ) derived an emp irical calibration fo r R23 b ased on Tg-metallicities for a sample of H II regions. This 
calibration has been updated bv iPilvugin fc ThuanI (|2005l . ; hereafter P05), using a larger sample of H II regions. 
They perform fits to the relationship between R23 and Tg-metallicities that includes an excitation parameter P that 
corrects for the effect of ionization parameter. The resulting calibration has an upper branch calibration that is valid 
for Tg-based metallicities 12 -I- log(0/H)> 8.25, and a lower branch calibration that is valid for Tg-based metallicities 
12 -I- log(0/H)< 8.0. We use the [N II]/[0 II] ratio (FigurejBj) to discriminate between the upper and lower branches 
for P05, and we apply the appropriate upper and lower-branch calibrations (equations 22 and 24 in P05): 
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-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 
log ([Nll]/Hct) log{([NII]/Ha)/([OIII]/HpJ} 

Fig. 10. — The observed relationship between the metallicities derived using the Tg method and (a) the [N II] A6584/Ha Une ratio, and 
(b) the [O III]/H/3/[N II]/Ha ratio. SDSS galaxies with useable [O III] A4363 line fluxes (S/N> 3) are shown as black filled circles. The 
PP04 calibration (dashed lines) and 95 percentile lines (dotted lines) are shown for each line ratio. The D02 calibration (red dot-dashed 
line) is shown for panel (a). A fraction of SDSS-Te galaxies have Te metallicities that lie below the 95 percentile line in both the [Nil] /Ha 
and [O III]/H/3/[N II]/Ha diagrams. These galaxies have high [NII]/Ha and [NII]/[OII] ratios, suggesting high (above solar) metallicities. 
The Te-method appears to underestimate the metallicity in these galaxies, possibly as a result of temperature gradients or fluctuations 
that may occur preferentially at high metallicities. 



19-ui ^n/m fl23 + 726.1 + 842.2P + 337.5P^ 

+ log(U/Hj„ppe, - g^^gg ^ ^^^^p ^ ^ l.793i?23 ^ ^ 



, , i?23 + 106.4+ 106. 8P- 3.40P2 

12 + bg 0/H io^er = ^ (A12 

8V / flower 17.72 + 6.60P + 6.95P2 - 0.302^23 

where P23 = ^^3727,29+[Om] AA4959,5007 ^ p ^ [OIII]AA4959.5007/H/3 ^^^.^^^^ ^^^^ ^^^^^^^^ 

reproducing Te-based metallicities with the P05 calibration is ~ 0.1 dex 

Because the P05 method was derived using fits to H II regions, we test whether the P05 method is applicable to the 
relationship between the SDSS Te metalhcities and R23 in Figure [TTl In Figure [TTb. we plot all SDSS galaxies in our 
sample with S/N ratio > 3 in the [O III] A4363 line. The upper and lower P05 branches are shown for different values 
of the excitation parameter P (red dot-dashed lines). Several galaxies lie outside the bounds of the P05 lower branch 
(12 -h log(0/H) (Te)< 8.0). 

In Figure [TTb . we exclude the galaxies that lie below the lowest 95 percentile line in the PP04 03N2 calibration 
(Figure \Wb) that are likely to have unreliable Te metallicities. As we discussed in Section IA.4.11 these excluded 
galaxies have [N II] /Ha and [N II] /[O II] ratios that indicate metallicities above solar. We note that the excluded 
galaxies have excitation parameters between 0.2 < P < 0.8, with a mean excitation parameter of 0.46 ± 0.03. These 
P values are more consistent with the P05 upper branch (range 0.2 < P < 0.8, mean 0.64 ± 0.03) than the P05 lower 
branch (range 0.6 < P < 1.0, mean 0.8 ± 0.1) for our SDSS sample. The Te method may not be reliable for these 
galaxies. 

Combined Te-strong-line method 
Demcolo, Terlevich & Terlevich (2002) - D02 

The iDenicolo et al.l (|2002f ) calibration is based on a fit to the relationship between the Te metallicities and the 
[N II] /Ha line ratio for ~ 155 H II regions. Of these H II regions, ^ 100 have metallicities derived using the Te 
method, and 55 H II have metallic i ties e stimated using the theoretical 1VI91 R23 method, or an empirical method 
proposed bv lDiaz &: Perez-lVIonterol (|2000 f) method based on the sulfur lines. The division between H II regions with 
Te-based metallicities and those with strong-line metallicities occurs at 12 -|-log(0/H)~ 8.4. The D02 calibration is 
given by a linear least-squares fit: 

12 + log(0/H) = 9.12 + 0.73 X Ar2 (A13) 

where N2 — log([NII] A6584/Ha). D02 estimate that the uncertainty the derived metallicities is ^ 0.2 dex. 

In Figure [TUl we compare the D02 fit (red dot-dashed line) to the [N II]/Ha-Terelationship for the SDSS galaxies. 
The D02 method provides a reasonable fit to the SDSS galaxies, given the large scatter, and is similar to the PP04 N2 
curve to within ~ 0.2 dex over the metallicity range 7.4 < 12 + log(0/H) < 8.8. 
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Fig. 11. — The observed relationship between the metaUicities derived using the Te method and the R23 Une ratio for (a) all SDSS 
galaxies in our sample with measurable (S/N> 3) [O III] A4363 fluxes, and (b) for the SDSS galaxies in our sample with measurable 
[O III] A4363 lines that lie above the lowest 95 percentile line in the PP04 03N2 calibration (Figure [TOj . 

METALLICITY CONVERSIONS: WORKED EXAMPLES 

Three galaxies have metaUicities of 12 + log(0/H)= 8.3,8.6, and 9.1 calculated using three different methods; 
KK04, PP04, and D0 2, respectively. To compare these galaxy metaUicities with those derived by the SDSS team 
(jTremonti et al.ir2004f ). we convert the three galaxy metaUicities into a metallicity base of T04. 

For a galaxy wit h meta llicity 12 + log(0/H)/f/fo4 = 8.3, the KK04 metallicity is calculated from the lower R23 
branch (see Section IA.2.3P . Table [3] gives the coefficients of the polynomial that converts a KK04 metallicity into a 
T04 base metallicity. Because our KK04 metallicity is from the lower branch, we use the linear relation for the lower 
branch conversion: 

[log (0/H) + 12]to4 = -4.5710 + 1.53261 x [log (0/H) + 12] 
-8.2 

To convert an original PP04 metaUicity (12 + Iog(0/H)ppo4 = 8.6) from 03N2 into a T04 base metallicity. Tabled 
gives 

[log (0/H) + 12]to4 = -738.1193 + 258.96730 x [log (0/H) + 12]ppo4 
-30.057050 X [log (0/H) + 12]|,po4 
+1.167937 X [log (0/H) + 12]|,po4 
-8.9 

An original D02 metallicity 12 + log(O/H)p)02 = 9.1 cannot be converted into a T04 metallicity because this D02 
abundance is above the valid range for our conversion from D02 into T04 (8.05 — 8.9). 

As a final example, we convert a T04 abundance of 12 + log(O/H)T04 ~ 8.3 into an 1VI91 base metallicity. Note that 
the valid upper and lower branch ranges overlap for the conversion of T04 into M91, ie. 8.2 < log(O/H)T04 < 9.2 
(upper branch) and 8.05 < log(O/H)T04 < 8.4 (lower branch). At a T04 metallicity between 8.2 < log(O/H)T04 < 8.3, 
the M91 R23 calibration is reaching a local maximum and is insensitive to metallicity. In this regime, the M91 upper or 
lower branch should be selected based on the [N II] /[O II] or [N II] /Ha ratio (Section lA.ip . If the log([N II] /Ha) ratio 
in our example galaxy is —0.7, then Figure [H indicates that the metallicity is on the R23 upper branch, and therefore, 
the conversion from a T04 metallicity of 12 + log(O/H)T04 = 8.3 into an M91 base metallicity is 



[log (0/H) + 12]m91 =404.1716- 131.53250 : 



log (0/H) 
2 



12] 



T04 



-14.49175 X [log (0/H) + 12]j,o4 



-0.5285842 x [log (0/H) 
'8.6 



12]|^o4 
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TABLE 1 

Comparison of the 10 metallicity calibrations 



# 


ID 


Emission lines 


Calibration Class 


Reference 


1 


T04'^ 


[Oil], H/3, [OIII], Ho, [Nil], [SII] 


theoretical 


Tremonti et al. C2004~) 


2 


Z94 


R23 


theoretical 


Zaritskv ct al. (1994) 


3 


KK04 


R23, [OIII] /[Oil] 


theoretical 


Kobulnickv & Kcwlci^ (2004') 


4 


KD02 


[NII]/[OII], R23, [OIII]/[OII] 


theoretical 


Kewlcv & Dopita (2002) 
McGauKh (19917 


5 


M91 


R23, [OIII] /[Oil] 


theoretical 


6 


D02 


[Nil] /Ha 


combined 


Dcnicolo ct al. (2002) 


7 


PP04 


[Nil] /Ha, [OIII]/H/3 


empirical 


Pcttini & Paecl (2004) 


8 


PP04 


[Nil] /Ha 


empirical 


Pcttini & Pagcl (2004) 


9 


P01,P05 


R23, [oni]/[oii] 


empirical 


Pilvudn (2001): Pilvudn & Thuan (2005) 


10 


Te 


[OIII] A4363, [OIII] AA4959, 5007 


direct 


AUer (1984): Stasihska (2005): Izotov et al. (2006a) 



^ The T04 method uses a statistical technique to calculate the probability distribution 
of an object having a particular metallicity based on model fits to the [Oil], H/3, [OIII], 
Ha, [Nil], [SII] emission lines. 



TABLE 2 

Robust fits'' to the MZ relations for the 9 strong-line metallicity 

calibrations 



ID a b c d rms residuals 



T04 


-0.694114 


1.30207 


0.00271531 


-0.00364112 


0.12 


Z94 


72.0142 


-20.6826 


2.22124 


-0.0783089 


0.13 


KK04 


27.7911 


-6.94493 


0.808097 


-0.0301508 


0.10 


KD02 


28.0974 


-7.23631 


0.850344 


-0.0318315 


0.10 


M91 


45.5323 


-12.2469 


1.32882 


-0.0471074 


0.11 


D02 


-8.71120 


4.15003 


-0.322156 


0.00818179 


0.08 


PP04 03N2 


32.1488 


-8.51258 


0.976384 


-0.0359763 


0.10 


PP04 N2 


23.9049 


-5.62784 


0.645142 


-0.0235065 


0.09 


POl 


91.6457 


-25.9355 


2.67172 


-0.0909689 


0.12 


P05 


41.9463 


-10.3253 


1.04371 


-0.0347747 


0.13 



^ Robust fits are of the form y — a bx -\- hx^ + bx^ where y — log(0/H) + 12 and 
X — log(M) where M is the stellar mass in solar units 



TABLE 3 

Metallicity Calibration Conversion Constants' 



ID KK04 Z94 M91 KD02 T04 D02 PP04 03N2 PP04 N2 

y (x) (x) (x) (x) (x) (x) (x) (x) 



KK04 

x-range'' 8.4-9.3 8.4-9.1; 8.0-8.25 8.4-9.2; 8.05-8.3 8.3-9.2 ; 8.05-8.4 8.2-8.9 ; 8.05-8.4 8.2-8.9 ; 8.05-8.3 8.2-8.9 ; 8.05-8.3 

branch'^ up; low up; low up; low up; low up; low up; low 

a 348.1710 -355.2968; 2.0021 1149.479; 1.0944 -162.1918; 5.0521 -1.7278; 6.5483 389.1568; 6.5339 211.1405; 5.2218 

6 -117.65370 114.8835; 0.77696 -389.9349; 0.87842 57.57935; 0.39887 1.52612; 0.21508 -133.79070; 0.21761 -81.11275; 0.37813 

c 13.502640 -12.09838; 44.33080; -6.533486; -0.034389; 15.59350; 10.577470; 

d -0.5130281 0.4258375; -1.675994; 0.249956; ; -0.6020232; -0.4510567; 

Pr 0.003 0.007; 0.007 0.031; 0.004 0.048; 0.029 0.065; 0.039 0.041; 0.041 0.064; 0.038 

Z94 

x-range'' 8.55-9.2 8.4^9.1 8.4-9.2 8.4^9.2 8.05-8.9 8.05-8.9 8.05-8.8 

a -1112.0910 -868.280 1086.903 11.2595 63.6386 230.9335 40.5515 

b 379.11370 291.6262 -366.5700 -1.47641 -13.87785 -76.73906 -8.67461 

c -42.880040 -32.42779 41.41521 0.136681 0.872216 8.711059 0.581691 

d 1.6218750 1.206416 -1.554630 -0.3244087 

Pr 0.005 0.010 0.040 0.056 0.081 0.052 0.076 

■M91 

x-iaiigc'= 8.25-9.15 8.4^9.2 8.4-9.2 ; 8.1-8.4 8.2-9.2 ; 8.05-8.4 8.2-8.9 ; 8.05-8.4 8.2-8.9; 8.05-8.4 8.2-8.8 ; 8.05-8.3 

branch'^ up; low up; low up; low 8.2-8.9; 8.05-8.4 up; low 

a -641.2458 393.9855 890.1.334: -1.0907 404.1716; 3.0825 85.2839; 4.7927 267.3936; 2.4196 87.3710; 0.7641 

6 226.42050 -127.86040 -295.90760; 1.12114 -131.53250; 0.61779 -18.63342; 0.40796 -85.20014; 0.70167 -19.39959; 0.90362 

c -26.374200 14.050330 33.004210; 14.491750; 1.130870; 9.219665; 1.193544; 

d 1.0268860 -0.5109532 -1.2227730; -0.5285842; ; -0.3267103; ; 

Pr 0.008 0.008 0.034; 0.004 0.047; 0.044 0.071; 0.059 0.048; 0.076 0.064; 0.070 

■Kno2 

x-range'= 8.2-9.2 8.4-9.2 8.5-9.1 ; 8.05-8.3 8.2-9.2 ; 8.05-8.4 8.2-8.9 ; 8.05-8.4 8.2-8.9 ; 8.05-8.3 8.2-8.9 ; 8.05-8.3 

branch'^ up: low up: low up; low up; low up; low 

a 75.5327 -476.8939 -2127.7470; 1.0007 387.2871: 4.1582 47.3054; 5.9875 159.0567; 5.8961 1094.5410:4.5323 

6 -23.64323 160.45270 720.35980; 0.88853 -129.69190; 0.49751 -10.06952; 0.27338 -54.18511; 0.28562 -388.67530; 0.45232 

c 2.658484 -17.752170 -81.051380; 14.718030; 0.649502; 6.395364; 46.233760; 

d -0.0948578 0.6579900 3.0435570; -0.5531547; ; -0.2471693; -1.8276270; 

Pr 0.038 0.036 0.038; 0.003 0.041; 0.039 0.048; 0.055 0.047; 0.059 0.047; 0.054 

T04 

x-range'' 8.6-9.15 ; 8.2-8.4 8.4^9.2 8.6-9.1 ; 8.05-8.4 8.1-9.2 8.05-8.9 8.05-8.9 8.05-8.9 

branch'^ up; low up; low 

a -461.2352; -4.5710 -472.8841 -69.7024; -1.0200 1.3782 193.9000 -738.1193 -1661.9380 

b 158.44840; 1.53261 158.20310 16.68313; 1.13063 0.52591 -64.87895 258.96730 585.17650 

c -17.946070; -17.414190 -0.880678; 0.036003 7.411102 -30.057050 -68.471750 

d 0.6828170; 0.6427315 ; -0.2756653 1.1679370 2.6766690 

Pr 0.060; 0.060 0.062 0.062; 0.063 0.051 0.072 0.058 0.072 

D02 

x-range'' 8.2-9.2 8.4-9.3 8.5-9.1 ; 8.0-8.4 8.05-9.2 8.05-9.2 8.05-8.9 8.05-8.9 

branch'^ up; low 

a 1202.5280 -114.3143 -121.3958; -1.0361 680.5167 253.0031 -1.3992 -629.0499 

6 -412.30820 32.79523 28.86410; 1.13468 -235.27350 -87.03697 -5.32702 215.37940 

c 47.332730 -2.782591 -1.599613; 27.345830 10.241880 1.562757 -24.305910 

d -1.8065160 0.0731175 : -1.0552390 -0.3984731 -0.0938063 0.9168766 

Pr 0.058 0.056 0.061; 0.093 0.038 0.044 0.033 0.000 

PP04 (03N2) 

x-range'' 8.2-9.2 8.4-9.3 8.5-9.1 ; 8.05-8.4 8.1-9.2 8.05-9.2 8.05-8.9 8.05-8.8 

branch'^ up: low 

a 631.2766 52.2389 -65.0991; 2.1063 664.8453 230.7820 36.6598 -8.0069 

6 -210.02090 -18.67559 15.74995; 0.74427 -225.75330 -75.79752 -7.64786 2.74353 

c 23.483050 2.447698 -0.837514; 25.768880 8.526986 0.508480 -0.093680 



TABLE 3 — Continued 



ID 


KK04 


Z94 


M91 


KD02 


T04 


D02 


PP04 03N2 


PP04 N2 


y 


(x) 


(x) 


(x) 


{x) 


(x) 


[x) 


(x) 


(x) 


0, 


-U.o 1 U^ZoD 


-U. lUl 10 1 O 




-U.y ( D±ODo 


-u.oiozoy'i 








pr 


U.UOU 


U.U4 ( 


U.UOD, U.U ( o 


U.U4 i 


U.U4D 


U.U40 




n riQfl 


PP04 (N2) 


















x-range'' 


8.2-9.2 


8.4-9.3 


8.5-9.1 ; 8.05-8.4 


8.05-9.2 


8.05-9.2 


8.05-8.9 


8.05-8.9 




branch'^ 






up; low 












a 


916.7484 


656.5128 


1334.9130; 3.1447 


569.4927 


319.7570 


-444.7831 


512.7575 




b 


-309.54480 


-224.11240 


-464.86390; 0.61788 


-192.51820 


-107.13160 


165.42600 


-180.47540 




c 


35.051680 


25.734220 


54.166750; 


21.918360 


12.208670 


-20.202000 


21.415880 




d 


-1.3188000 


-0.9812624 


-2.0986640; 


-0.8278840 


-0.4606539 


0.8249386 


-0.8427312 




Pr 


0.068 


0.065 


0.071; 0.050 


0.045 


0.051 


8.557e-6 


0.042 





^ The conversion constants convert metallicities from the calibrations given in row 1 (x) into a metalUcity that is consistent with the calibration given in column 1 (y), using y — a-\-bx-\-cx^ -\-dx^ . The 
range over which a metallicity using a calibration given in row 1 (a:) can be converted into a metalUcity that is consistent with the calibration given in cokimn 1 (y).^ Some conversions are different, 
depending on whether the input metallicity {x) is on the upper or lower branch (up; low). In some regimes, both the upper and lower branch conversion equations are valid. In these cases, 
additional information such an [N II]/ [O II] or [N II] /Ha ratio is needed to distinguish between the upper and lower R23 branches. A worked example for this case is given in AppcndixfB] 



TABLE 4 
Residual Metallicity Discrepancy 



Base 


Mean Residual 


ID* 


Discrepancy'' (dex) 


KK04 


0.011 


Z94 


0.018 


M91 


0.011 


KD02 


0.014 


T04 


0.034 


D02 


0.034 


PP04 03N2 


0.012 


PP04 N2 


0.021 



a b 

calibration into which the other 7 MZ relations have been converted Mean residual metalUcity discrepancy between the 7 converted MZ 

relations and the base MZ relation. 



TABLE 5 

Relative median metallicity and rms scatter, for, the 9 str,onc;-line 
metallicity calibrations from 30,000 random sets of sdss cjalaxies 



ID T04 Z94 KK04 KD02 M91 D02 PP04 (03N2) PP04 (N2) POl Mean 

med (rms) med (rms) med (rms) med (rms) med (rms) med (rms) med (rms) med (rms) med (rms) rms 



T04 •■■ 0.006 (0.08) 0.053 (0.09) 0.033 (0.07) 0.036 (0.09) 0.076 (0.10) 0.040 (0.08) 0.053 (0.09) 0.046 (0.14) 0.09 

Z94 -0.006 (0.08) • ■ ■ 0.046 (0.05) 0.023 (0.05) 0.024 (0.04) 0.062 (0.09) 0.029 (0.06) 0.042 (0.08) 0.034 (0.11) 0.07 

KK04 -0.053 (0.09) -0.046 (0.05) • • • -0.023 (0.05) -0.014 (0.02) 0.016 (0.06) -0.017 (0.05) -0.003 (0.07) -0.004 (0.09) 0.06 

KD02 -0.033 (0.07) -0.023 (0.05) 0.023 (0.05) ••• 0.005 (0.05) 0.039 (0.06) 0.006 (0.05) 0.019 (0.05) 0.011 (0.10) 0.06 

M91 -0.036 (0.09) -0.024 (0.04) 0.014 (0.02) -0.005 (0.05) ••• 0.031 (0.08) -0.001 (0.06) 0.009 (0.07) 0.010 (0.08) 0.06 

D02 -0.076 (0.10) -0.062 (0.09) -0.016 (0.06) -0.039 (0.06) -0.031 (0.08) ••• -0.033 (0.05) -0.020 (0.03) -0.017 (0.11) 0.07 

PP04 (03N2) -0.040 (0.08) -0.029 (0.06) 0.017 (0.05) -0.006 (0.05) 0.001 (0.06) 0.033 (0.05) ••• 0.014(0.05) 0.012(0.12) 0.07 

PP04 (N2) -0.053 (0.09) -0.042 (0.08) 0.003 (0.07) -0.019 (0.05) -0.009 (0.07) 0.020 (0.03) -0.014 (0.05) ••• 0.004 (0.11) 0.07 

POl -0.046 (0.14) -0.034 (0.11) 0.004 (0.09) -0.011 (0.10) -0.010 (0.08) 0.017 (0.11) -0.012 (0.12) -0.004 (0.11) ••• 0.11 



